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1.1 INTRODUCTION 

There are many differences and also remarkable similarities between the ion 
chemistry and physics of planetary ionospheres and the ion chemistry and 
physics of astronomical environments beyond the solar system. 

In the early Universe, an expanded cooling gas of hydrogen and helium was 
embedded in the cosmic background radiation field and ionized by it. As the 
Universe cooled by adiabatic expansion, recombination occurred and molecular 
formation was driven by catalytic reactions involving the relict electrons and 
protons. Similar chemical processes are effective in the ionized zones of gaseous 
and planetary nebulae and in stellar winds where the ionization is due to 
radiation from the central stars, in the envelopes of supernovae where the 
ionization is initiated by the deposition of y-rays, in dissociative shocks where 
the ionization arises from electron impacts in a hot gas and in quasar broad-line 
region clouds where the quasar is responsible for the ionization. At high 
altitudes in the atmospheres of the Jovian planets, the main constituents are 
hydrogen and helium and the ion chemistry and physics is determined by the 
same processes, the source of the ionization being solar ultraviolet radiation 
and cosmic rays. 

After the collapse of the first distinct astronomical entities to emerge from 
the uniform flow, heavy elements were created by nuclear burning in the cores 
of the collapsed objects and distributed throughout the Universe by winds and 
explosions. The chemistry and physics became more complicated. Over 90 
distinct molecular species have been identified in interstellar clouds where they 
are ionized globally by cosmic ray impacts and locally by radiation and shocks 
associated with star formation and evolution. Complex molecules have also 
been found in circumstellar shells of evolved stars. At intermediate and low 
altitudes in the Jovian atmospheres, the ion chemistry is complicated by the 
increasing abundance of heavy elements such as carbon, and an extensive array 
of complex molecules has been predicted. Reactions involving heavy elements 
dominate the structure of the ionospheres of the terrestrial planets and the 
satellites Titan and Triton. 
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1.2 HYDROGEN AND HELIUM PLASMAS 

We begin with a discussion of the chemistry and physics of various astrophysical 
and planetary plasmas in which the major gas components are hydrogen and 
helium. After establishing the elements of the hydrogen-helium plasma chem- 
istry, we will consider the complexities that arise with the introduction of heavier 
elements. 


1.2.1 ASTROPHYSICAL ENVIRONMENTS 

Much of the Universe consists of plasmas dominated by hydrogen and helium. 
We compare the chemistry and physics of several such plasmas. 


1.2. 1.1 Early Universe 

The formation of the first molecules in the Universe was a result of chemical 
reactions in which electrons and protons acted as catalysts. According to the 
standard cosmology, space occupied by the Universe expanded adiabatically. 
At t = \0 4 s, the temperature was 10 12 K. The Universe then cooled rapidly 
as it expanded, reaching a temperature of 10 9 K after 100 s. At 10 9 K, protons 
and neutrons had formed and a brief period of nucleosynthesis occurred during 
which a substantial abundance of helium w r as produced together with trace 
amounts of deuterium and ^Li. The matter was embedded in a radiation field 
characterized by the same temperature as the matter. The radiation field 
contained energetic photons and photoionization ensured that all the atoms 
were fully stripped. 

As the expansion continued, the supply of energetic photons diminished as 
the temperature continued to decrease. After some tens of thousands of years, 


Li 3+ and He 2 + captured electrons by radiative recombination 

Li 3+ + e * Li 2 + + hv (1) 

He 2+ -he ► He + + hv (2) 

while the reverse processes of photoionization 

hv + Li 2+ ► Li 3+ 4- e (3) 

hv + He + ► He 2+ + e (4) 

were ineffective. The initial recombinations were soon followed by 

Li 2+ + e ► Li + + hv (5) 

He + -h e ► He + hv (6) 
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and Li + and He could survive photoionization by the radiation field because 
of their high ionization potentials of respectively 75.62 and 24.56 eV. 

With the formation of neutral helium, the first molecular ions He 2 + and 
LiHe + came into existence through the processes of radiative association 

He + + He ►He 2 + +/«v (7) 

Li + + He >LiHe + +/iv (8) 

In the radiative association of He + and He, during their approach along 
the 2 I g potential energy curve, a photon may be emitted in an electronic 
transition to bound vibrational levels of the ground 2 £ u + state of He 2 . The 
rate coefficient of reaction (7) has been calculated (Stancil et al. 1993). The 
radiative association of Li + and He is a much less probable event because the 
radiative transition involved is a spontaneous emission from the vibrational 
continuum of the ground ‘I state into the discrete vibrational levels. 

The molecular ion He 2 + has a dissociation energy of 2.47 eV ( Ackermann 
and Hogreve, 1991) and an effective threshold for ground state photodissocia- 
tion near 115 nm. The vibrational level populations are in thermal equilibrium 
at the cosmic background temperature and the ions are quickly removed by 
photodissociation (Gaur et al, 1988) 

He 2 + + hv ► He + + He (9) 

by the background radiation field. They are also removed by dissociative 
recombination 

He 2 + + e ►He + He (10) 

and by charge transfer 

He 2 + + H ► He + He + H + (11) 

The rate coefficients of reaction (10) and (11) are uncertain. 

The ion LiHe + is removed by photodissociation into the excited electronic 
‘I and 'n states 

LiHe + + hv >Li + +He (12) 

and by dissociative recombination 

LiHe + + e > Li + He (13) 

The rate coefficient of reaction (13) is unknown. It may be slow compared to 
most such processes. In any event, photodestruction severely limits the abun- 
dances of He 2 + and LiHe + that can be achieved and the ions play an 
insignificant role in the early chemistry. 

The presence of neutral helium leads also to the formation of HeH + by 
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radiative association 


H + + He * HeH + + hv (14) 

in a process analogous to reaction (8). The rate coefficient is small with a value 
at 10000 K of 1.3 x 10 20 cm 3 s“ 1 (Roberge and Dalgarno, 1992). Once 
formed, HeH + is only slowly destroyed. Photodissociation by absorption into 
the vibrational continuum of the ground state is negligible (Roberge and 
Dalgarno, 1992) and photodissociation by absorption into the excited A 1 ! state 
has a high threshold near 19eV. Dissociative recombination 

HeH + + e ► He + H (15) 

which proceeds with a rate coefficient measured by Yousif and Mitchell (1989) 
to be 1 x 10“ 8 (300/T) 1/2 cm 3 s“\ and 

HeH + + H ► H 2 + + He (16) 

which proceeds with a rate coefficient of (9.1 ± 2.5)10“ 10 cm 3 s~ 1 (Karpas et 
aL , 1979), are the main destruction mechanisms. 

The processes involving helium and lithium ions were preliminary events to 
the main activity which began when the temperature had decreased to about 
4000 K and photoionization of atomic hydrogen 

H + /iv ► H + + e (17) 

had ceased to be effective. Then hydrogen was rapidly neutralized by radiative 
recombination 

H + + e ► H -f hv (18) 

The Universe entered the recombination era which was also the dawn of 
chemistry, giving rise to the formation of the first neutral molecule H 2 . The 
particle density at this time of 10 5 years was about 10 3 cm' 3 . 

Direct formation of molecular hydrogen by radiative association of ground 
state hydrogen atoms is highly forbidden because it involves a triplet-singlet 
transition from the b 3 E u + state of H 2 to the X 1 I 9 + state. The suggestion has 
often been made that the excited state population may be sufficient such that 

H(2s) + H ► H 2 + hv (19) 

H(2p) + H ► H 2 + hv (20) 

which involve allowed transitions, make significant contributions to the forma- 
tion of H 2 . Latter and Black (1991) have made estimates of the rate coefficients. 
The process (20) can be viewed alternatively as a Raman scattering of Lyman 
alpha radiation from the vibrational continuum of the ground electronic state, 
terminating in vibrational levels of the ground state. Raman scattering, as a 
mechanism for forming H 2 , was investigated by Federman and Frommhold 
(1982). 
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More important is the process of associative ionization 

H(n = 2) + H ► H 2 + + e (21) 

Urbain et al. (1991) have presented the cross sections for H(2s). The possible 
influence of associative ionization has not yet been included in any model of 
the early Universe, but it has been incorporated in the envelope chemistry of 
Supernova 1987a (Miller et al, 1992; Yan et al, 1994). 

With hydrogen not entirely neutral, H 2 + molecular ions could be formed by 
radiative association 

H + + H ► H 2 + + hv (22) 

and negative hydrogen ions H " could be formed by radiative attachment 

H+e * H~ + hv (23) 

The reverse reactions of photodissociation 

H 2 + + hv ► H + + H (24) 

and photodetachment 

H ~+hv >H+e (25) 


delayed but could not prevent the formation of hydrogen molecules by charge 
transfer 


H 2 + + H — 

-*H 2 + H + 

(26) 

and associative detachment 



H + H- 

— ► H 2 + e 

(27) 

In these sequences, the protons and 
formation of H 2 . 

Mutual neutralization processes 

the electrons act 

as catalysts for the 

H 2 + + H“ — 

-►H + H + H 

(28) 

H + + H" — 

— ► H + H 

(29) 

He + + H — 

— ► He + H 

(30) 

play a minor role in limiting the abundance of H 2 . 

Molecular hydrogen can be destroyed by (Chupka 
Maujean, 1989; Dalgarno, 1993a) 

et al., 1968; Glass- 

H(n = 2) + H 2 

► H 3 + + e 

131) 

It is reformed in part by dissociative recombination 


Hj + + e — 

~^H 2 + H 

(32) 

but 

H 3 + + e 

+ H + H 

(33) 
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may also be a significant channel. The rate coefficient of dissociative recombina- 
tion of H 3 + is a crucial parameter in chemical models of interstellar clouds, 
supernova envelopes and planetary ionospheres. It has experienced a checkered 
history. The most recent experiments (Yousif et al. , 1991; Canosa et al ., 
1992; Larsson et al. , 1993) appear to establish a rate coefficient of order 
10 _7 (300/T) 1 2 cm 3 s 1 but the actual mechanism responsible for the fast rate 
has yet to be established (Bates, 1992; Bates et al. , 1993). 

An additional source of H 2 + is reaction (16), but the reverse process 

H 2 + + He ► HeH + 4- H (34) 

which is exothermic for vibrational levels of H 2 + with v > 3 also occurs. 

Some chemistry took place involving the trace constituents deuterium and 
lithium. Lithium remained mostly ionized because of the scarcity of electrons 
left over after hydrogen recombination. Some formation of LiH + occurred 
through (Lepp and Dalgarno, 1987) 

Li + + H ► LiH + + hv (35) 

and, more important, formation of HD through 

D + + H 2 ► H + + HD (36) 

The reaction (36) is the major source of HD in interstellar clouds (Dalgarno et 
al ., 1973a). Its rate coefficient was measured by Ferguson (1986). 

Figure 1 shows the results of a calculation by Latter (1989) of the fractional 
abundances of the various forms of hydrogen and helium as a function of the 
redshift z. The relict ionization fraction after recombination is 1.8 x 10“ 4 . 



Figure LI The fractional abundance of the various ions of hydrogen and helium as a 
fraction of red shift. (From Latter, 1989) 
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1.2. 1.2 Gaseous nebulae and stellar winds 

In gaseous nebulae and planetary nebulae, a warm ionized zone is created by 
the action of the central stars. The main components of the gas are hydrogen 
and helium. With increasing distance from the ionizing sources, geometrical 
dilution and optical depth effects reduce the ionizing flux and the gas becomes 
largely neutral. Some ionization persists from the photoionization of minor 
constituents C, S and Si. The formation of molecular hydrogen occurs through 
the reaction sequences (22) to (27) (Black, 1978). Of greater interest is the 
formation of the molecular ion HeH + because it may be detectable through its 
infrared vibrational emission. In gaseous nebulae, its formation is enhanced by 
the radiative association He + and H 

He + + H ♦ HeH + + hv (37) 

which has a rate coefficient of about 10“ 15 cm 3 s“ 1 (Zygelman and Dalgarno, 
1990). Reaction (37) is not effective in the early Universe because of the small 
overlap of ionized helium and neutral hydrogen. Figure 1.2 is an illustration 
of the distributions of H“, H 2 + , H 2 and HeH + in a nebula of density 
n(H) = 10 4 cm" 3 and n(He) = 10 3 cm “ 3 as functions of distance from a central 
star with a hydrogen ionizing flux of 10 48 s“‘ and a blackbody temperature of 



Figure 1.2 The abundances of H 2 , H“, H 2 + and HeH + in a model spherical nebula 
with a density oflO 4 cm 3 corresponding to a central star with an effective temperature 
of 200000 K as a function of radial distance. (From Cecchi-Pestillini and Dalgarno, 1993) 
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200000 K (Cecchi-Pestillini and Dalgarno, 1993). A search for HeH + in the 
planetary nebula NGC 7027 (Moorhead et ai , 1988) was unsuccessful. Model 
calculations (Cecchi-Pestillini and Dalgarno, 1993) suggest that the emission 
may be close to the observational upper limit. 

High-velocity winds are generated by newly formed protostars. The winds 
may be partly neutral (Natta et ai , 1988) and molecular synthesis may be driven 
by reactions of the ionized and neutral components in a chemistry closely 
analogous to that occurring in the early Universe (Rawlings et ai , 1988; 
Glassgold, 1991). The gas is heated by ultraviolet photons from the protostar. 
In addition to the H + and H " sequences, three-body recombination 

H + H + H *H 2 + H (38) 

contributes to the formation of molecular hydrogen in the higher density winds. 
Reaction (38) was important also in the evolution of the first collapsing objects 
in which cooling due to molecular hydrogen controlled the fragmentation mass 
spectrum (Stahler et ai , 1986). 

1.2.1. 3 Supernova 1987a 

The envelope of Supernova 1987a is composed mostly of hydrogen and helium 
with a minor component of heavy elements appropriate to the interstellar 
medium of the Large Magellanic Cloud galaxy. The gas is heated and ionized 
by events initiated by the absorption of y-rays emitted in the decay of 60 Co. 
The ion chemistry is modified by an important additional sequence for the 
formation of H 2 in which the radiative association of He + and H in reaction 
(37) is followed by reactions (16) and (26). 

The He + ions can be removed by reaction with H 2 . It appears that despite 
the substantial exoergicity, He + reacts very slowly with H 2 in the lowest v = 0 
and v = 1 vibrational levels but very rapidly for H 2 in vibrational levels greater 
than or equal to v ~ 2 (Jones et ai , 1986). The corresponding rate coefficient 
is 1 x 10 9 exp(— 12000/T) cm 3 s~ The reaction is probably (Jones et ai , 
1986) 

He + + H 2 (r > 2) ► He + H + H + (39) 

Reactions (21) and (31) involving excited hydrogen atoms are significant sources 
of H 2 + and H 3 + in the envelope. 

The chemistry of the supernova envelope is of special interest because of the 
possible identification of emission features at 3.41 and 3.53 pm as transitions 
between rotation-vibrational levels of H 3 + (Miller et ai , 1992). Figure 1.3 
illustrates the match obtained between the observational data (Meikle et ai , 
1989) and a thermal emission spectrum at a temperature of 2050 K and an 
expansion velocity of 3200 km s~ 1 (Miller et ai , 1992). The identification was 
stimulated by the positive identification of H 3 + emission in the spectra of Jovian 
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Day 192 (Meikle et aJ.,1989) 
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Figure 1.3 The points are observational data on the infrared emission spectrum of 
supernova 1987a. The solid line is a theoretical curve showing the thermal emission 
spectrum at 2050 K and an expansion velocity of 3200 km s' 1 for H 3 + . (From Miller 
et ai ., 1992) 


auroras (Trafton et ai , 1989; Drossart et ai , 1989, 1993; Kim et ai , 1990; 
Maillard et ai , 1990; Miller et ai , 1990; Oka and Getalle, 1990; Geballe et ai , 
1993; Trafton et ai , 1993). 

1.2. 1.4 Quasar broad-line regions 

Broad emission lines in allowed transitions observed toward quasars appear to 
originate in small dense high-velocity clouds. They are ionized by the ultraviolet 
and X-radiation from the quasar, and the formation of molecular hydrogen in 
the broad-line region clouds proceeds by the same reactions that occur in the 
supernova (Kallman et ai , 1987; Crosas and Weisheit, 1993). 

1.2. 1.5 Dissociative shocks 

Fast shocks result from protostellar outflows, stellar winds and supernova 
explosions in the interstellar gas. If the shock velocity exceeds 50 km s' 1 , all 
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log N h { cm' 2 ) 


Figure 1.4 The abundance of various ions of hydrogen and helium in gas of density 
10 5 cm~ 3 subjected to a shock of velocity 80 km s~ 1 as a function of column density of 
hydrogen behind the shock. (From Neufeld and Dalgarno, 1989) 


the molecules are destroyed in the heated gas and ionization is produced by 
electron impacts. The ionized gas then recombines and cools and, as it does, it 
releases ultraviolet photons. This precursor radiation modifies the chemistry 
behind and ahead of the shock (Hollenbach and McKee, 1989; Neufeld and 
Dalgarno, 1989). The cooling gas eventually recombines and molecular hydro- 
gen is formed by reactions (22) and (26) and (23) and (27). The formation of 
molecular hydrogen occurs whilst the gas is still warm and it initiates chemistry 
in which the hydroxyl radical occupies a central role. If dust grains survive the 
shock, the abundance of molecular hydrogen may be enhanced by formation 
on the grain surfaces. Figure 1.4 illustrates the density profiles of the shocked 
gas as a function of the column density of shocked material (Neufeld and 
Dalgarno, 1989). 
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1.2.2 PLANETARY ENVIRONMENTS 

Most of the ionization in a planetary atmosphere occurs in a region called the 
thermosphere. The thermosphere is so named for the temperature structure of 
the terrestrial upper atmosphere, which increases rapidly from a local minimum 
(the mesopause) of 160-190 K near 85 km to a large and constant value (the 
exospheric temperature) of 700-2000 K, depending upon solar activity. This 
increase in temperature is caused by absorption of solar extreme and far 
ultraviolet radiation, while the tendency to approach a constant value results 
from the increase in conductivity that accompanies the exponentially decreasing 
number densities. In the terrestrial atmosphere, the base of the thermosphere 
occurs at a pressure of a few pbar. In the atmospheres of other planets, similar 
thermal structure is encountered, but the limiting temperature is not always 
large. Exospheric temperatures of 180-200 K have been reported for low solar 
activity on Mars (Nier and McElroy, 1977). On the nightside of Venus, the 
temperature actually decreases above the mesopause, from about 160 K to an 
exospheric temperature of about 100 K (Seiflf et ai , 1985), 

The same ultraviolet radiation that leads to thermospheric heating also 
produces ionization, and in this process energetic photoelectrons are ejected 
that cause further ionization. 

The chemical structure of the upper atmosphere is determined to a great 
extent by the location of the homopause (or turbopause), which is usually found 
in the lower thermosphere of a planet. Below the homopause, the major 
constituents are assumed to be thoroughly mixed by convection and/or turbu- 
lence. Above the homopause, diffusion becomes the dominant transport mech- 
anism, and the constituents separate out under the force of gravity. The mixing 
ratios of heavier species therefore decrease and those of lighter constituents 
increase with altitude. This simple distribution does not apply to minor species 
whose densities are photochemically determined, such as atomic oxygen in the 
terrestrial lower thermosphere. 

1. 2.2.1 Outer planets 

Unlike the terrestrial planets, the outer planets have retained their primordial 
atmospheres, the composition of which is similar to that of the sun. The 
dominant atmospheric species is H 2 , with a few percent He by number. There 
are also small quantities of methane and higher hydrocarbons present in the 
thermospheres of the major planets, but, because hydrocarbons are heavy, their 
densities decrease rapidly above the homopause. Photodissociation of H 2 and 
diffusive separation lead to increasing densities of hydrogen atoms above the 
homopause. At sufficient heights, atomic hydrogen is predicted to become the 
dominant constituent. Model atmospheres for the Jovian and Saturnian ther- 
mospheres are shown in Figures 1.5 and 1.6 (Smith et ai , 1983a; Kim, 1991). 
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Temperature (K) 

200 700 




Figure 1.5 (a) A model of the Jovian thermosphere. The altitude scale refers to height 

above the ammonia cloud tops (near 0.6 bar). The H 2 density and temperature profiles 
were taken from Atreya et ai (1981). The horizontal bars indicate the densities of H 2 , 
CH 4 and C 2 H 6 inferred from Voyager stellar occultation data. The other density profiles 
are computed without photochemistry, (b) Computed density profiles of hydrocarbons 
and hydrogen atoms. The values shown were computed with coupled hydrocarbon and 
ion chemistry. (Adapted from Kim, 1991) 




14 


A. DALGARNO AND J. L. FOX 


Temperature (K) 


100 200 300 400 500 



Figure 1.6 Model for the thermosphere of Saturn derived from the Voyager 2 solar 
and stellar occultation data. The altitude is referred to the 1 bar level. It should be noted 
that from the same data Festou and Atreya (1982) derived an exospheric temperature 
about twice that shown here. (Adapted from Smith et al . , 1983) 


They are consistent with the neutral density structure of these atmospheres 
inferred from the Voyager ultraviolet stellar occultation experiment (Smith et 
al , 1983a; Kim, 1991). 

Altitude profiles of the production rate of ions in the Jovian ionosphere are 
shown in Figure 1.7(a). Above 500 km, photoionization of H 2 

H 2 + *v * H 2 + + e (40) 

is the major ion production process. Below 500 km, photoelectron impact 
ionization 

e + H 2 ► e + H 2 + + e (41) 

assumes a greater importance, because of the steepening of the incident 
spectrum as the radiation penetrates into the atmosphere. This radiation 
hardening leads to energetic photoelectrons. Although H 2 is the major atmo- 
spheric constituent between 350 and 450 km, photoionization of atomic hydro- 
gen is the most efficient source of ionization because the photons that are 




ION CHEMISTRY IN ATMOSPHERIC AND ASTROPHYSICAL PLASMAS 15 




Figure 1.7 (a) Computed production rate profiles of ions in the Jovian ionosphere for 

solar maximum conditions and a solar zenith angle of 22.5 . The hydrocarbon ions are 
labeled in order of increasing peak production rate. The curve labeled e-H 2 + is the 
production rate of H 2 + from electron impact ionization of H 2 . (b) Computed ion density 
profiles for a local time of noon for a latitude of 22.5 C. A magnetic dip angle of 45 ° 
is assumed. Altitude is referred to the ammonia cloud tops. (Adapted from Kim, 1991) 
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energetic enough to ionize molecular hydrogen are absorbed at higher altitudes. 
Photons with wavelengths between the ionization threshold of H 2 (i; = 0) at 
80.4 nm and the ionization threshold of H at 91.2 nm penetrate down and ionize 
H in the photoionization process (17). 

Although H 2 + is the major ion produced over a large altitude range of the 
Jovian ionosphere, it is quickly transformed to H 3 + in the fast ion-molecule 
reaction 

H 2 + -F H 2 >H 3 + 4- H (42) 

In steady state, the density of H 2 + is accordingly very small. 

Protons are produced in direct ionization of hydrogen atoms and in dissocia- 
tive ionization of hydrogen molecules. Radiative recombination (18) is slow 
with a rate coefficient of 4 x 10“ 12 (250/7^)°‘ 7 cm 3 s" \ where T e is the electron 
temperature (Bates and Dalgarno, 1962). The protons may diffuse downwards 
and be removed at high densities by the three-body recombination (Miller et 
aU 1968 ) 

H + + H 2 + H 2 >H 3 + + H 2 (43) 

producing H 3 + . The charge transfer of H* to H 2 in the ground vibrational 
state is endothermic by 1.8 eV, but it becomes exothermic for vibrational levels 
v > 4. Then (McElroy, 1973a) 

H + + H 2 (v > 4) ► H + H 2 + (44) 

may be an efficient mechanism for removing ionization that could help in 
resolving the discrepancy between the electron densities inferred from the radio 
occultation experiments on the Pioneer 10 and 11 and Voyager 1 and 2 
spacecraft and the theoretical models (Chen, 1981; McConnell et al 1982; 
Waite et al , 1983; Atreya, 1986; Cravens, 1987; Majeed et al , 1991). The 
measured electron density peaks are an order of magnitude higher and occur 
at an altitude 600 km above the theoretical predictions made neglecting vi- 
brationally excited H 2 . 

The electron densities in Figure 1.7(b) were calculated assuming a vibrational 
temperature of 1500 K over the entire altitude range. Charge transfer to 
vibrationally excited H 2 is not effective under all conditions and ion transport 
along field lines driven by meridional winds or electric fields has been invoked 
(Atreya, 1986; Majeed and McConnell, 1991). 

Similar problems have been encountered in models of the ionosphere of 
Saturn. A possible resolution will be discussed in Section 1.3.2. 

The process (44) may contribute to the removal of protons in photon- 
dominated regions of the interstellar medium, but it is there counter-balanced 
by the photoionization of vibrationally excited H 2 (Black and van Dishoeck, 
1987) 


H 2 (y > 4) + hv 


H 2 + + e 


(45) 



ION CHEMISTRY IN ATMOSPHERIC AND ASTROPHYSICAL PLASMAS 17 


The H 3 + molecule is produced by reactions (42) and (43) and removed 
at high altitudes by the dissociative recombination reactions (32) and (33). 
Emission from H 3 + in the infrared has been observed in auroras on Jupiter 
(Trafton et ai y 1989; Drossart et al 1989, 1993; Kim et al 1990; Maillard et 
al , 1990; Miller et al. , 1990; Oka and Geballe, 1990), Uranus (Trafton et al , 
1993) and Saturn (Geballe et al. , 1993). Figure 1.8 is an illustration of the 
spectrum and the fit that can be obtained by assuming thermal equilibrium at 
a temperature of HOOK (Miller et al. , 1990). The excitation processes have 
been discussed by Kim et al. (1992). 

There is a significant production of Fie + ions in the high ionospheres of the 
outer planets. The ions recombine radiatively as in reaction (6) but slowly. 
Radiative association (37) and radiative charge transfer (Jura and Dalgarno, 
1971; Zygelman and Dalgarno, 1990) 

He + + H ► He + H + + hv (46) 

neutralize the helium ions but the reaction (39) with vibrationally excited 
molecular hydrogen may be a more efficient He + removal process. 

In models of the Jovian ionosphere it has been assumed that the rate 
coefficient of the reaction of He 4 with H 2 is equal to the upper limit 
2 x 10“ 13 cm 3 s 1 obtained in room temperature measurements (Boringer and 


Hg -1- emission spectrum around 2.1 um. 
Temperature = 1100 K; /p=0.56 



Figure 1.8 The points are the observational data on the infrared emission spectrum of 
Jupiter. The solid curve is the thermal emission spectrum at 1 100 K predicted for H 3 + . 
(From Miller et al 1990) 



18 


A. DALGARNO AND J. L FOX 


Arnold, 1986) and that, in addition to the end products in reaction (39), channels 
leading to H 2 + and HeH + also occur (McElroy, 1973a; Johnsen el al , 1980; 
Atreya, 1986). If reaction (39) is adopted instead, the predicted He + concentra- 
tions are larger, by as much as two or three orders of magnitudes. 

At altitudes near the homopause, the abundance of H + and He + are greatly 
modified by the presence of hydrocarbon molecules in the atmosphere. 


1.3 PLASMAS WITH AN ADMIXTURE OF 
HEAVY ELEMENTS 

The introduction of even a small fraction of heavy elements into a hydro- 
gen-helium plasma drastically modifies its chemical composition, ionization 
balance and thermal structure. In interstellar clouds, carbon, oxygen and 
nitrogen make up at most one part in a thousand by number of the material, 
yet the millimeter and radiofrequency spectrum is dominated by molecules into 
which one or more heavy elements has been incorporated. In the atmospheres 
of the outer planets below the altitudes of the homopause, the increased relative 
abundances of heavy atoms leads to a diverse array of complex molecular 
species. Methyl acetylene has been detected on Jupiter and possibly on Saturn. 

1.3.1 ASTROPHYSICAL ENVIRONMENTS 

Other than the early Universe, astrophysical plasmas contain an important 
component of heavy elements and their chemistry involves an elaborate set of 
ion-molecule reactions. 

1.3.1. 1 Diffuse and translucent interstellar clouds 

An interstellar cloud is a collection of material in interstellar space in which 
the constituent particles move with the same mean velocity. Diffuse and 
translucent molecular clouds are identified by their absorption of the light from 
individual stars. Absorptions by the diatomic molecules CH, CN, CH + , OH, 
CO, H 2 , C 2 and NH have been observed in the visible, ultraviolet and infrared. 
The diatomic molecule CS has been found through its emission in the millimeter 
region. The absence of polyatomic systems is ascribed to the destructive effects 
of the interstellar radiation field. The temperatures of diffuse clouds, at least of 
their central cores, is low, of the order of 30 K. Because the clouds are cold, 
the chemistry is driven by ionization. Translucent clouds are more extensive 
versions of diffuse clouds through which optical and infrared radiation can still 
penetrate. 

Ionization is created by interstellar ultraviolet photons, which produce C + , 
S + , Si + and ionized metals, and by cosmic rays which ionize everything. The 
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carbon chemistry initiated by ultraviolet ionization starts with the radiative 
association 


C + + H 2 >CH 2 + +/iv (47) 

followed by dissociative recombination 


CH 2 + + e — 

— ♦ CH + H 

(48) 

to produce CH. The CH 2 + ion may 
abstraction reaction 

react with H 2 to form CH 3 + 

in the 

ch 2 " + h 2 - 

— » CH 3 + + H 

(49) 

The dissociative recombination 



CH 3 + +e- 

— ► CH + H 2 

(50a) 

— 

— ♦ CH, + H 

(50b) 

leads to CH and CH 2 . In dense clouds 
is small, radiative association 

in which the fractional electron density 

av + h 2 — 

— CH 5 + + hv 

(51) 

may occur before dissociative recombination. 

The CH 2 , made in reaction (50b), will undergo photodissociation 


CH 2 + hv- 

— » CH + H 

(52) 

to again produce CH. 

In diffuse clouds, CH is removed by photodissociation 


CH + hv - 

— ► C + H 

(53) 

and by photoionization 



CH + hv — 

-*CH f +e 

(54) 

It is also converted to CO and CN by the neutral exchange reactions 


CH + O — 

-» CO + H 

(55) 

CH + N — 

— ► CN + H 

(56) 

The bringing together of two carbon atoms to form C 2 can be accomplished 
by the sequence 

C + CH — 

— >C, 4 + H 

(57) 

C 2 * + H,— 

— ♦ C 2 H + + H 

(58) 

C 2 H + + e- 

— C 2 + H 

(59) 

c 2 + + c — 

-*C 2 + C' 

(60) 
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Figure 1.9 The carbon chemistry in interstellar clouds 


The carbon chemistry is summarized in Figure 1.9. 

Branching reactions with oxygen and nitrogen, (55), (56) and 

C 2 + 0 >C O + C (61) 

C 2 + N ► CN + C (62) 

produce CO and CN. 

The sequence 

CH 3 + + O > HCO + + H 2 (63) 

HCO + + e ► H + CO (64) 

is also a source of CO. A similar sequence with N leads to CN. The CN 

molecules can be converted to CO by 


CN + O » CO + N (65) 

Carbon monoxide is the end product of numerous sequences, because of its 
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large dissociation energy. It is lost by photodissociation 

CO T hv ► C T O (66) 

and by the reaction with He + 

He + T CO * He + C + + O (67) 

Because the reaction of He + with H 2 is slow at low temperatures, reaction (67) 
is an efficient destruction mechanism in cold clouds. In warmer regions, reaction 
(39) removes the helium ions. 

This system of reactions rests upon the radiative association of C + with H 2 . 
Models of interstellar clouds (Black and Dalgarno, 1977; Black et al , 1978; 
van Dishoeck and Black, 1986, 1989; Viala, 1986; Viala et al , 1988) require a 
rate coefficient of about 5 x 10“ 16 cm 3 s -1 for the process. An experimental 
upper limit of 1 x 10“ 15 cm 3 s -1 has been obtained (Luine and Dunn, 1989; 
Gerlich and Hornung, 1992). 


Cosmic ray ionization is needed in diffuse clouds to produce 0 + 
The 0 + ions are produced mostly by charge transfer of H + 

and N + . 

H + + O — 

— * H + 0 + 

(68) 

They react with H 2 in a series of abstraction reactions 


o + + h 2 - 

— ♦ OH + + H 

(69) 

OH + + H,- 

— ♦ H 2 0 + + H 

(70) 

h 2 o + + h 2 - 

— * H 3 0 + + H 

(71) 

terminating in H 3 0 + . Then dissociative recombination 


H,0 + + e- 

— * H 2 0 + H 

(72a) 

- 

— * OH + H 2 

(72b) 


leads to H 2 0 and OH. The branching ratios for reaction (72b) have been 
measured recently (Adams, 1992) in what is a major advance for the interstellar 
chemistry of dense clouds. The branching ratio is less critical in diffuse cloud 
chemistry because any H 2 0 formed is rapidly photodissociated 

H 2 0 + hv * OH + H (73) 

to yield OH. The sequence is illustrated in Figure 1.10. 

The formation of OH provides other pathways for production of CO 


C + + OH — 

— CO + H + 

(74) 

C + + OH — 

— CO + + H 

(75) 

CO + + h 2 — 

— HCO + + H 

(76) 

HCO + + e — 

— H + CO 

(77) 
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Figure 1.10 The oxygen chemistry in interstellar clouds 


Cosmic ray ionization is particularly effective in producing OH because it 
does so through the ionization of the major constituents H and H 2 to produce 
H + and H 2 + . The H 2 + ions are quickly converted to H 3 + by reaction (42) 

H 2 + 4- H 2 ► H 3 + + H 


(78) 
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and H + and H 3 + then react with O atoms. The analogous reactions do not 
occur with N atoms and the nitrogen ion molecular chemistry is initiated by 
the direct cosmic ray ionization of N to produce N + . The reaction 

N + + H 2 ► NH + + H (79) 

is slightly endothermic (Gerlich and Hornung, 1992) and in cold clouds must 
compete with 

N + + CO * NO + C + (80) 

► N 4- CO + (81) 

In warmer diffuse clouds, reaction (79) proceeds and is followed by an abstrac- 
tion sequence leading to NH 3 . Photodissociation of NH 3 and NH 2 leads to 
NH. Photodissociation of NH itself 

NH + hv ► N + H (82) 

is fast (Kirby and Goldfield, 1991) and it is difficult to reproduce the observed 
abundance of NH by gas phase chemistry (Crutcher and Watson, 1976; Meyer 
and Roth, 1991). Formation of NH 3 , NH 2 and NH on the surfaces of grains 

and their release into the gas phase provides a possible resolution of the 

discrepancies (Wagenblast et al , 1993). 

The detection of CS in a diffuse cloud (Drdla et al. , 1989) is significant in 
providing information on the sulfur chemistry. The ionization structure of 
molecular clouds and perhaps the sizes of collapsing regions of gas are sensitive 
to the chemistry of sulfur (Hartquist et al ., 1993). Like carbon, sulfur is ionized 
by the interstellar radiation field 

S + hv S + + e (83) 

Charge transfer from C + 

C + +S >C + S + (84) 

is an additional source of S\ Like C + , S + does not react with hydrogen but 
it can undergo radiative association 

S + + H 2 * H 2 S + + hv (85) 

The rate coefficient of reaction (85) is uncertain. A value of 10 _17 (T/300) 0 1 
cm 3 s -1 (Miliar and Herbst, 1990) is usually adopted. Because the process is 
slow, ionization due to S + persists deeply into interstellar clouds. The slowness 
of reaction (85) permits the efficient formation of CS through (Drdla et al 1989) 


S + + CH — 

— CS + + H 

(86) 

s * + c 2 — 

— cs + + c 

(87) 

cs + + h 2 — 

-> HCS + + H 

(88) 

HCS + + e — 

— > H + CS 

(89) 
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Photodissociation of CS is the main destruction channel. The exothermic 
reaction 


CS + O ► CO + S (90) 

has an activation energy of 760 K (Millar et al, 1991) and is slow in diffuse 
clouds. It is effective in warm gas. Other reactions that can destroy CS tend to 
enter sequences that lead back to it. 

The chemical models that have been constructed of diffuse clouds, and indeed 
of translucent and dense clouds, are based on sets of reactions and rate 
coefficients listed by Anicich and Huntress (1986), Millar et al (1991) and 
Anicich (1993). For reactions for which no measurements exist it is assumed 
that reactions of ions with homonuclear molecules occur at the Langevin rate 
of about 10 9 cm 3 s \ independent of temperature. For reactions of ions with 
heteronuclear molecules, one or another of several theoretical models (Clary, 
1990; Dubernet et al , 1992; Wickham et al., 1992) is adopted. In agreement 
with experiments (Adams et al ., 1985; Marquette et al ., 1985; Clary, 1990), the 
theoretical models indicate a rate coefficient increasing with decreasing tem- 
perature. The theoretical models also demonstrate that the rate coefficients are 
sensitive to the rotational level of the molecule. In astronomical environments, 
the rotational populations are often far from thermal equilibrium. Reaction rate 
coefficients may also be affected by the fine-structure populations in the case 
of ions such as N + (Marquette et al , 1988) and C + and Si + (Graff, 1989). 

The models have met with mixed success (Brown and Rice, 1986; Millar et 
al , 1987; Suzuki et al , 1988; Herbst and Leung, 1989, 1990) despite the 
flexibility offered by the uncertainties in the values of the reaction rate 
coefficients at the ambient cloud temperatures. Perhaps large molecules partici- 
pate in the chemistry. There are observational evidence and theoretical argu- 
ments (Leger and Puget, 1984; Allamandola et al , 1985; Lepp et al, 1988) to 
support the view that a substantial component of polycyclic aromatic hydro- 
carbons and fullerenes is present in the diffuse interstellar medium. In diffuse 
clouds the large molecules will be singly and doubly ionized (d’Hendecourt and 
Leger, 1987; Lepp and Dalgarno, 1988a; Millar, 1991), at least in the boundary 
region, and interesting chemistry may occur (Lepp and Dalgarno, 1 988b; Millar, 
1991). 

Translucent clouds are intermediate in extent between diffuse and dense 
clouds (van Dishoeck and Black, 1980; Lepp and Dalgarno, 1988b; van 
Dishoeck, 1990). Similar to translucent clouds is a class of high-latitude clouds, 
seen through CO millimeter emission (van Dishoeck and Black, 1980). In some 
translucent and high-latitude clouds, the polyatomic molecules H 2 CO and 
C 3 H 2 have been detected in addition to the diatomic molecules observed in 
diffuse clouds. Formaldehyde can be made from (Dalgarno et al, 1973b) 


CH 3 + 0 


H 2 CO + H 


(91) 
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and perhaps 

CH 5 + + O * H 3 CO + + H 2 (92) 

H 3 CO + + e * H 2 CO -f H (93) 

but not in the abundances that are found. Similarly cyclopropenylidene can be 
made by a sequence of reactions involving C + , but again not in the observed 
abundances. Large molecules may play a role. The vulnerability of H 2 CO and 
C 3 H 2 to photodissociation means that an efficient production mechanism is 
needed. 

None of the models is successful in reproducing the abundances of CH + , 
which is readily destroyed by dissociative recombination 


CH + + e- 

— > C + H 

(94) 

photodissociation 

CH + + hv — 

— ♦ C + H + 

(95) 

and reactions with H and H 2 

CH + + H — 

-c + + h 2 

(96) 

CH + + H 2 — 

- ch 2 + + H 

(97) 


A major additional source of CH f is needed. It seems necessary to invoke 

C + + H 2 > CH+ + H (98) 

Reaction (98) is endothermic by 0.4 eV so it will proceed rapidly only if the 

ambient gas is warm or the H 2 vibrational levels are anomalously excited or 
the C + ion is accelerated (Duley et a /., 1992). It must also occur in a gas with 
a high H/H 2 ratio to avoid the overproduction of CH that would otherwise 
result. Shocked gas appears to be excluded by the absence of a velocity 
difference between CH + and CH or CO. 

However the CH + is produced, it may be an important source of CO and 
CN through reactions following 

CH + H- O ► CO + 4- H (99) 

CH + + N > CN + + H (100) 

Because the models underestimate CH + , they may also underestimate CO and 
CN. 

Many diffuse molecular clouds may be the boundary regions of more 
extensive clouds. When subjected to intense radiation fields they are called 
photon-dominated regions or photodissociation regions (PDRs). They occur in 
the nuclei of external galaxies and starburst galaxies (Stacey et ui , 1991) and 
in the Milky Way, the Orion Molecular Cloud being a dramatic example (Stacey 
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et ai , 1993). The chemistry is driven by photoionization and limited by 
photodissociation (Tielens and Hollenbach, 1985a, 1985b; Sternberg and 
Dalgarno, 1989). The irradiated edge of the PDR may be heated to temperatures 
above 1000 K and the chemistry is modified. Photodissociation can be more 
efficient than dissociative recombination in removing molecular ions. 


1.3. 1.2 Dense molecular clouds 

Dense molecular clouds are the birthplaces of stars. Emission from molecules 
formed within them cools the clouds and brings regions of them to the point 
of gravitational collapse. The cooling radiation is detected and provides a 
diagnostic probe of star formation. Table 1.1 (Dalgarno, 1993b) is a collection 
of molecular species that have been detected in interstellar clouds. Although 
ion-molecule chemistry in itself does not explain the observed abundances, 
ion-molecule reactions play a major part in determining the chemical composi- 
tion of dense clouds. 

The ion-molecule chemistry begins with the production of H 2 + by cosmic 
rays followed by reaction (42). Proton transfer reactions such as (69) followed 
by abstraction sequences (68) to (70) lead to polyatomic molecular ions which 
undergo dissociative recombination to form neutral molecules. Figure 1.10 
illustrates the reaction sequence. 

Mechanisms for building the complex hydrocarbons of Table 1.1 are illu- 
strated in Figure 1.11. The radiative association process 

CH 3 + + H 2 >CH 5 + + hv (101) 


is a critical step. There remains a considerable uncertainty in its rate coefficient 
at the temperatures of dense clouds near 10-20 K (Gerlich and Hornung, 1992). 
Insertion reactions such as 


c + + ch 4 — 

-c 2 h 3 + + h 

(102) 

followed by 



C 2 H 3 + + e — 

— C 2 H 2 + H 

(103) 

makes acetylene, which reacts again with C + 


c + + c 2 h 2 — 

— C 3 H + + H 

(104) 

to continue the build-up of carbon atoms. Condensation reactions 


ch 3 + + ch 4 - 

— ♦ c 2 h 5 + + h 2 

(105) 


also help to form complex organic molecules. 

Rate coefficients for radiative association are generally uncertain. Yet the 
processes are an essential component of the formation of the complex organic 
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Table 1.1 

Interstellar molecules 



H 2 

Hydrogen 

CH 

Methylidyne 

CH + 

Methylidyne ion 

OH 

Hydroxyl 

C, 

Carbon 

CN 

Cyanogen 

CO 

Carbon monoxide 

NO 

Nitric oxide 

cs 

Carbon monosulfide 

SiO 

Silicon monoxide 

so 

Sulfur monoxide 

NS 

Nitrogen sulfide 

SiS 

Silicon sulfide 

PN 

Phosphorous nitride 

HC1 

Hydrogen chloride 

SiN 

Silicon nitride 

H,0 

Water 

C 2 H 

Ethynyl 

HCN 

Hydrogen cyanide 

HNC 

Hydrogen isocyanide 

HCO 

Formyl 

HCO + 

Formyl ion 

N 2 H" 

Protonated nitrogen 

H,S 

Hydrogen sulfide 

HNO 

Nitroxyl 

OCS 

Carbonyl sulfide 

SO, 

Sulfur dioxide 

HCS + 

Thioformyl ion 

c,o 

Carbon suboxide 

C 2 S 

Dicarbon sulfide 

h 2 co 

Formaldehyde 

h 2 cs 

Thioformaldehyde 

NH 3 

Ammonia 

HCNS 

Isothiocyanic acid 

HNCO 

Isocyanic acid 

HOCO + 

Protonated carbon dioxide 

C 3 H 

Propynylidyne 

c 3 n 

Cyanoethynyl 

c 3 s 

Tricarbon sulfide 

c 3 o 

Tricarbon monoxide 

c 2 h 2 

Acetylene 

H,0 + 

Hydronium ion 

HCNH + 

Protonated hydrogen cyanide 

C 3 H 2 

Cyclopropenylidene 

ch 4 

Methane 

H.CCC 

Propadienylidene 

HCOOH 

Formic acid 

CH,CO 

Ketene 

hc 3 n 

Cyanoacetylene 

HNCCC 


HCCNC 

Ethynyl isocyanide 

c 4 h 

Butadinyl 

NH 2 CH 

Cyanamidc 

CH,CN 

Cyanomethyl radical 

CH ,NH 

Methanimine 

CH,CN 

Methyl cyanide 

H.CCCC 

Butatrienylidene 

CH,SH 

Methyl mercaptan 

C,H 

Pentynylidyne 

HCC 2 HO 

Propynal 

CH 3 OH 

Methyl alcohol 

ch 3 c 2 h 

Methyl acetylene 

NHXHO 

Formamide 

HC S N 

Cyanodiacetylene 

CH.CHCN Vinyl cyanide 

CH,NH 2 

Methylamine 

C h H 

Hexatrinyl 

hcooch 3 

Methyl formate 

CH ,CHO 

Acetaldehyde 

ch 3 c 4 h 

Methyl diacetylene 

ch',c,n 

Methyl cyanoacetylene 

ch,ch 2 cn 

Ethyl cyanide 

(CH,) 2 0 

Dimethyl ether 

ch 3 ch,oh 

Ethyl alcohol 

HC.N 

Cyanohexatriync 

HCqN 

Cyano-octatetrayne 

HC„N 

Cyanodecapentayne 




molecules found in dense clouds. Methyl alcohol may be a result of 


CH 3 + +H 2 0 * CH 3 + -H 2 0 + hv (106) 

CH , + ■ H 2 0 + e ►CHjOH + H (107) 

and ethyl alcohol a result of 

H 3 O l +C 2 H 4 ► H 3 0 + -C 2 H 4 + hv (108) 

HjO 1 • C 2 H 4 + c ► CH,CH 2 OH + H (109) 
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Figure 1.11 Mechanisms for the production of hydrocarbons in dense interstellar 
clouds 
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The hydrocarbons and their extensions illustrated in Figure 1.11 undergo 
branching reactions and incorporate other species. Reactions with N lead to 
the polyacetylenes. The reactions with O and with S are mostly unknown. The 
high relative abundances of C 2 S and C 3 S (Saito et al . , 1987; Yamamoto et al . , 
1987) suggest that sulfur insertion reactions 


S + C,H, + * HC 3 S + + H, 

(HO) 

are rapid (Smith et al , 1988; Suzuki et a /., 1988). The reactions with oxygen 
may be destructive (Bettens and Brown, 1992). They are often neglected. 

The neutral molecules are destroyed by reactions with He + . Acetylene, which 
is an important molecule in the hydrocarbon chain of Figure 1.11, is lost by 

He + + C 2 H, * He + H, + C, f 

(111a) 

♦ He + C 2 H 2 + 

(11 lb) 

» He + C 2 H + + H 

(111c) 

> He 4- CH + + CH 

(Hid) 

reaction (111b) being the major channel. Methane is also destroyed by He + 
in the reactions 

He + + CH 4 ♦ He + CH 4 + 

(112a) 

* He + CH, + + H 

(112b) 

» He + CH, + + H 2 

(114c) 

♦ He + CH + + H 2 + H 

(11 2d) 

* He + H + + CH, 

(1 12e) 


reaction (112b) being the major channel (Anicich, 1993). These reactions 
occur also in the ionospheres of the outer planets to which we now turn. 


1.3.2 OUTER PLANETS 

As shown in Figure 1.7(a), at an altitude below 350 km on Jupiter, a sharp layer 
of hydrocarbon ionization appears due to photons with wavelengths longer 
than the ionization thresholds of H, FL and He. The ionization potential of 
methane is 12.5 eV. Most larger hydrocarbons have still lower potentials. CH 3 
and can be ionized by the strong solar Lyman alpha line at 121.6 nm, 

corresponding to an energy of 10.2 eV. Photons are absorbed at the centers of 
the Lyman and Werner bands of H 2 but penetrate down in the wings, causing 
direct ionization of hydrocarbon molecules down to about 300 km (Kim and 
Fox, 1992). 
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Any He + ions are quickly removed by reactions such as (1 1 1) and (112). Of 
the H + ions that have diffused down, some can form H 3 + in the three-body 
reaction (43) but most react rapidly with CH 4 (Dheandhanoo et al . , 1984) 

H + + CH 4 ► CH 3 + + H 2 (113a) 

>CH 4 + + H (113b) 

and with C 2 H 6 (Mackay et al., 1981). 

On Saturn, the electron densities may be reduced to the measured values by 
postulating an influx of water (Chen, 1983; Shimizu, 1983; Connerny and Waite, 
1984; Majeed and McConnell, 1991). Protons can then be removed by 


H + + H 2 0 ► H + H 2 0 + 

(114) 

The H 2 0 + dissociatively recombines 

H 2 0 + + e ►OH + H 

or reacts with H 2 as in reaction (70), producing H 3 0 + . 

(115) 

The H 2 0 is also photodissociated 

H 2 0 + hv ► OH + H 

(116) 

and H + can charge transfer to OH 

H + + OH ► H + OH + 

Then reaction (70) and dissociative recombination 

(117) 

OH + + e ► O + H 

(118) 


occur. 

A detailed model of the ionosphere on Saturn (Majeed and McConnell, 1991) 
suggests that the presence of vibrationally excited H 2 or the injection of H 2 0 
from the rings at a flux exceeding 10 7 cm' 2 s _1 can reproduce the peak 
electron density, provided an upward drift along the field lines is imposed on 
the ions. In Figures 1.12 and 1.13, we present for the two theoretical models 
comparisons of the predicted ion and electron density profiles with the measure- 
ments recorded by the Voyager 2 spacecraft on the exit portion of the Saturn 
flyby. Little difference is found between the two models and the model topside 
scale heights are too large by a factor of three. In regions where transport is 
unimportant the scale height is inversely proportional to mass and directly 
proportional to temperature. Because it is unlikely that the ionosphere is colder 
by a factor of three or that the major ion is H 3 f and not H + , it appears it is 
not chemistry but dynamical effects due to strong winds or electric fields that 
control the ionization structure. 

In the region of the ionosphere near and below the Jovian homopause, direct 
photoionization of hydrocarbons and, to a lesser extent, charge transfer from 
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Figure 1.12 Computed profiles for H + and H 3 f and total electron densities in the 
Saturn ionosphere compared to Voyager 2 egress profile (curve a). The curves labeled 
H 3 + and H + were computed with an effective rate coefficient for loss of by reaction 
with vibrationally excited H 2 (reaction 44) of 1.38 x 10" 14 cm 3 s' 1 , while curve b was 
computed without this reaction. A vertical drift along field lines of 15 ms" 1 was imposed. 
(Adapted from Majeed and McConnell, 1991) 


H + and He + to neutral hydrocarbon molecules produce a layer of hydrocarbon 
ions. Figure 1.14 shows that the major source of hydrocarbon ions with a single 
carbon atom below the homopause is photoionization of CH 4 and CH 3 , and 
the major loss processes near 320 km are reactions with methane and acetylene, 
such as (105) and 

CH 5 + + C 2 H 2 >C 2 H 3 + +CH 4 (119) 



Figure 1.13 Computed profiles for H + and H 3 + and total electron densities in the 
Saturn ionosphere compared to Voyager 2 egress profile (curve a). The curves labeled 
H 3 + and H 1 were computed with an influx of water from the rings of 2.2 x 10 7 
cm" 2 s' 1 , while curve b was computed without such a flux. A vertical drift along field 
lines of 15 ms' 1 was imposed. (Adapted from Majeed and McConnell, 1991) 
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Figure 1.14 Sources (a) and sinks (b) of Cj hydrocarbon ions in the Jovian ionosphere, 
(a) The labels pi(CH 4 + ) and pi(CH 3 + ) refer to the photoionization of CH 4 and CH 3 
leading to CH 4 + and CH 3 + , respectively. The labels 1, 2 and 3 are the production rates 

due to the reactions H 3 + + CH 4 ► CH 5 + + H 2 , H + + CH 4 ► CH 4 + + H and 

H + + CH 4 ► CH 3 + + H 2 , respectively, (b) The labels rec(CH 3 + ) and rec(CH 5 + ) 

refer to dissociative recombination of CH 3 + and CH 5 + , respectively. The labels 1 and 

2 refer to the reactions CH 5 + +C 2 H 2 — ^C 2 H 3 + + CH 4 and CH 3 + + CH 4 > 

C 2 H 5 + ■+■ H 2 . The dotted curve is the total rate of reactions that convert into C 2 
hydrocarbon ions. (Adapted from Kim, 1991) 


The largest sources of hydrocarbon ions with two carbon atoms near 320 km 
are photoionization of C 2 H 2 , C 2 H 4 and C 2 H 6 and the most rapid loss processes 

C 2 FV +c 2 h 2 — >c 3 h 3 + + ch 4 
c 4 h 5 + + h 2 


are 


c 2 h 5 + + c 2 h 2 


( 120 ) 

( 121 ) 
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and to a lesser extent reactions with CH 4 such as 

C 2 H 5 + + CH 4 ► C 3 H 7 + + H 2 (122) 

Above about 380 km, most of the hydrocarbon ions are destroyed in dissociative 
recombination but near 320 km almost all the hydrocarbon ions produced react 
to form C 3 and C 4 hydrocarbon ions. Rate coefficients for possible reactions 
of large hydrocarbon ions are not well known, but it is likely that the large 
ions will also react, especially with acetylene, to produce higher hydrocarbon 
ions which will ultimately be destroyed by dissociative recombination. 

The major ions in the ionospheres of the outer planets and the reactions that 
connect them are shown in Figure 1.15. In one model of the Jovian E region 
(Kim and Fox, 1992), reactions of hydrocarbons with one or two carbon atoms 
were included explicitly. Higher hydrocarbons were taken into account as tw'o 
groups of ions called C 3 H n + and C 4 H n + , which were assumed to be lost by 
dissociative recombination with rate coefficients of 3.5 x 10" 7 (300/Te)° 5 
cm 3 s~\ Figure 1.16 shows computed altitude profiles of electron densities and 
the major ions H f , H 3 \ He + and H 2 + , and the two groups of hydrocarbon 
ions for a latitude of 22.5 ° and a local time of noon. The peak density of 
hydrocarbon ions is greater than 10 4 citU 3 near 320 km. 

Although the solar flux at Uranus is 400 times less than at Earth, electron 
density profiles measured by the radio occultation experiment on Voyager 2 
show' a number of layers peaking below about 2000 km (Lindal et ai , 1987). 
An influx of H 2 0 or CH 4 from the rings could reduce the electron densities in 
the equatorial regions (Strobel et ai , 1991) in much the same way as predicted 
for Saturn (Connerny and Waite, 1984). Because Uranus does not have as strong 
an internal heat source as do Jupiter and Saturn, the tropopause is very cold, 
about 50 K. Consequently, the hydrocarbon densities at higher altitudes are 
much smaller (Strobel et al ., 1991) and the hydrocarbon ion layer predicted for 
Jupiter and Saturn is probably absent. 

Electron density profiles determined from the Pioneer 10 and 1 1 radio 
occultation measurement of the Jovian ionosphere showed a number of layers 
below the main ionospheric peak (Fjeldbo et ai , 1974, 1976). The lower layers 
may be composed of metallic ions (Atreya et al. , 1974; Atreya and Donahue, 
1976), similar to those found in the terrestrial ionospheric E region (Atreya et 
ai , 1974), and which result from photoionization and charge exchange reactions 
of the major ions with metal atoms that are produced by ablation of meteorites 
(Hunten et ai , 1980; Hunten, 1981). The major meteoritic ions are Fe + , Mg + 
and Si + (Swider, 1992). The flux of meteors into the Jovian ionosphere may be 
enhanced over the terrestrial value by gravitational focusing (Grebowsky, 1981). 
For plausible meteoroid fluxes at Jupiter in contrast to the terrestrial case, 
photoionization is not important compared to charge exchange as a source 
of Fe + ions and the peak metallic ion production rates are similar to those of 
H 2 + , H + , He + and CH 3 + . The peak densities were estimated to be in the 
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Figure 1.15 Chemical network among ions in the Jovian ionosphere. The major 
production processes for the species shown are indicated by solid arrows and the major 
loss processes are shown by dashed arrows. (Adapted from Kim, 1991) 
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Figure 1.16 Computed ion density profiles for H + , H 3 + , He + and H 2 T and two groups 
of hydrocarbon ions: C,H n + and C 4 H„ 4 for a local time of noon and a latitude of 
22.5 (Adapted from Kim, 1991) 


10 3 -10 6 cm 3 range, with similar contributions from other metal ions, such as 
Mg + . The observed layers, which have thicknesses of about 100 km, can be 
reproduced by postulating fluxes of Na + and S + ions from the Galilean 
satellites and high wind shears (Chen, 1981). The models described here, 
however, considered only radiative recombination as a loss mechanism for 
metal ions. As in the terrestrial ionosphere, the ions may form clusters which 
are destroyed by dissociative recombination. 


1.4 HEAVY ELEMENT PLASMAS 

By heavy element plasma we mean plasmas in which hydrogen and helium play 
a lesser role, though they may still have a significant impact on the behavior 
of the plasma. The ionospheres of the terrestrial planets and the satellites Titan 
and Triton are examples, as is the mantle of the Supernova 1987a. 

1.4.1 SUPERNOVA 1987a 

The detection of carbon monoxide in the ejecta of Supernova 1987a about 100 
days after the explosion of a blue supergiant star in the Large Magellanic Cloud 
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testifies to the durability of molecular systems in environments that are 
apparently inimical to their formation and hostile to their survival. 

The ejecta are subjected to y-rays and fast positrons produced in the 
radioactive decay of 5b Ni to 56 Co to 56 Fe. At later times, the longer-lived 
isotopes 57 Co and 44 Ti become the main sources of continuing energy 

deposition into the ejecta. The y-rays undergo Compton scattering with free 
and bound electrons and are degraded into X-rays. The X-rays are preferentially 
absorbed by heavy elements and by Auger ionization produce multiply charged 
positive ions. The positive ions capture electrons in radiative 

X m+ -he ► x (m " 1)+ -h hv (123) 

and dielectronic 

X" + -he > X <W1 + 1) + * (124a) 

► X (ffl - 1)+ + hv (124b) 

recombination. 

Once some neutral material has been produced, charge transfer processes 
such as 


X m+ + H ♦ X (m ' n+ 4 H + (125) 

X m+ -h He * X (w “ 1)+ 4 He + (126) 

quickly reduce the ionization stages and the plasma becomes a mixture of 
electrons, singly charged ions and neutral atoms. A tentative identification of 
an emission feature at 22.92 pm to Fe 2+ has been made (Moseley et ai, 1989). 
The presence of Fe 2+ would not be expected because of the rapidity of the 
charge transfer recombination processes (Dalgarno, 1990) 

Fe 2+ + He * Fe + 4 He + (127) 

Fe 2+ 4 0 > Fe + 4 0 + (128) 

The distribution of the ionization after the charge transfer of multiply charged 
ions has occurred depends on the abundances of the elements and the rate 
coefficients for the charge transfer of singly charged ions in collision with neutral 
atoms. The values are uncertain. Of particularly importance are the reactions 
of He + ions. Calculations of the potential energy curves of the HeO + (Augustin 
et ai y 1973) indicate that the most probable mechanism for the charge transfer 

He + 4 O ♦ He 4 0 + (129) 

is spin-orbit and rotational coupling of 2 n and 2 I molecular states. A rate 
coefficient of the order of 10“ 13 cm 3 s“ 1 is possible. It cannot be much greater 
because of the observation that He + ions exist in abundance in the topside 
ionosphere. 
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A value of 3.3 x 10 9 cm 3 s 1 has been recommended for 

He + + Si * He + Si + (130) 

but it appears to be an arbitrary estimate (Prasad and Huntress, 1980). 

Explicit calculations of the potential energy curves and coupling matrix 
elements of HeC + have been carried out. The charge transfer 

He + + C ► He + C + (131) 

appears to be driven by a rotational coupling between the D 2 Z” and B 2 n 
states, leading to a rate coefficient at 3000 K of 1 x 10" 13 cm 3 s" 1 (Kimura et 
ai y 1993). 

In contrast to the reactions of He + , charge transfer of ions of complex atoms 
are likely to be fast because of the accessibility of favorably located potential 
energy curves. If so, the ionization will flow from elements of high ionization 
potential like oxygen to elements of low ionization potential like the metals Fe, 
Ni and Co in reactions such as 

0 + + C ► O + C + (132) 

C + + Fe ► C 4- Fe + (133) 

Estimates of the rate coefficients of charge transfer reactions between heavy 
ions and heavy atoms have been made (Pequignot and Aldrovani, 1986). 

The detailed distribution of the neutral atom and ion abundances depend on 
the specific charge transfer rate coefficients, the element abundances, the 
ultraviolet photon flux and the energetic electron flux, but the spectra are 
qualitatively consistent with a scenario dominated by exothermic charge trans- 
fer processes. 

At early times the formation of molecules is due to three-body processes such 
as 

He + C + O ► He + CO (134) 

but the gas is hot and the molecular species are removed by thermal dissociation 
processes such as 

He + CO ► He + C + O (135) 

e + CO ► e + C + O (136) 

Once the gas has cooled below about 5000 K, CO can survive and emission 
from vibrationally excited carbon monoxide was detected at 110 days (Oliva 
et ai , 1987; Rank et ai ., 1988; Spyromlio et al. , 1988). The 1-0 fundamental 
emission of SiO was also detected (Roche et al. , 1991; Wooden et al. , 1993). 
The diminishing intensity of the CO emission with the disappearance of 
emission from high vibrational levels could be attributed to the decreasing 
temperature. 
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Several mechanisms have been suggested that led to the formation of CO. 
At the time when CO was first detected, the particle density was about 
10 11 cm" 3 , too small for three-body collision to be effective. The initiating 
reactions had to be two-body collisions in which the formation of a composite 
system is stabilized by the emission of a photon. Radiative association 
(Petuchowski et al. , 1989; Lepp et al , 1990) 

C + -f O ► CO + + hv (137) 

followed by charge transfer 

CO + H-O >C0 + 0 + (138) 

is a possibility. The radiative association is slow with a rate coefficient of about 
10 18 cm 3 s” 1 (Dalgarno et al , 1990) and much of the CO + that is produced 
is lost by dissociative recombination 

CO + -he ► C + O (139) 

with a rate coefficient of 1 x 10" 7 (300/ T) 0 48 cm 3 s“ 1 (Mitchell and Hus, 1987). 
More efficient because the carbon is mostly neutral is the direct association 

C-hO ► CO -h hv (140) 

which has a larger rate coefficient of about 10” 17 cm 3 s -1 (Dalgarno et a /., 
1990). 

A potentially significant source is provided by sequences initiated by the 
formation of negative ions (Lepp et al , 1990). Thus O” is formed by radiative 
attachment 


O + e *0~+hv (141) 

with a rate coefficient of 10” 15 cm 3 s” L Then associative detachment 

C + O” * CO 4- e (142) 

leads to carbon monoxide. The efficiency of the sequence is impacted by the 
destruction of O by mutual neutralization 

He + + O” ► He + O (143) 

and photodetachment 

O - + hv ► O + e (144) 

Less direct sequences can make some contribution. The radiative association 

O + O ► 0 2 + hv (145) 

C + o 2 


followed by 


CO + o 


(146) 
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is an example. Still less direct but potentially significant in the core if the light 
elements are mixed back and surely important in the envelope is the sequence 
of reactions (37), (16) and (26) to form H 2 followed by 

C + H 2 ► CH + H (147) 

and reaction (55) or 

0 + H 2 ► OH + H (148) 

and 

OH + C ► CO + H (149) 

The radicals CH and OH can also be formed by radiative association. 

Satisfactory agreement with observations can be obtained, but it is necessary 
to limit the degree of microscopic mixing of helium with the carbon and oxygen 
(Liu et ai, 1992). 

The molecule SiO is also detected in the supernova ejecta (Roche et ai , 1991 ; 
Wooden et ai , 1993). It is made mostly by 

Si 4- O * SiO + hv (150) 

with some contribution from 

Si + 0 2 * SiO + O (151) 

It is destroyed by 

C + 4- SiO > Si + 4* CO (152) 

The rate coefficient of reaction (150) is unknown. A value of order 10“ 18 
cm 3 s“ 1 is sufficient to reproduce the derived abundances provided there is no 
fine-scale mixing of helium (Liu and Dalgarno, 1994). 


1.4.2 TERRESTRIAL PLANETS 
1. 4.2.1 Dayside ionospheres 

The atmospheres of Mars and Venus are composed mostly of C0 2 with small 
amounts of N 2 : 3.5% for Venus (von Zahn and Moroz, 1985) and about 2.7% 
for Mars (Owen et ai , 1977). The atmosphere of Mars also contains argon with 
an abundance comparable to that of the Earth, about 1.6%. In the thermo- 
spheres of both planets, photodissociation of C0 2 and N 2 results in the 
formation of small amounts of the photolysis products CO, O and N, along 
with photochemically produced 0 2 and NO. The compositions of the lower 
thermospheres of Venus and Mars are compared in Table 1.2 (Fox, 1992a). The 
major difference between the thermospheres of the two planets is the smaller 
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Table 1.2 Mixing ratios of minor constituents 
in the atmospheres of Venus and Mars (from Fox, 
1992b) 


Species 

Mars (120 km) 

Venus (135 km) 

O 

0.6( — 2)— 1 .2( — 2) 

7.3(-2) 

CO 

4.2( — 3) 

7.0< — 2) 

N 2 

2.5( — 2) 

7.4< — 2) 

NO 

7( — 5) 

3< — 5) 

o 2 

1 2.2( — 3) 

3( — 3) 

Ar 

1.5( — 2) 

6.3(-5) 

N 

2.1(— 4) 

1.6<-3) 

C 

2.5( — 5) 

4( — 4) 

He 


l(-6H2( — 6) 


mixing ratios of O and CO on Mars. Figures 1.17 and 1.18 show model 
thermospheres of Mars and Venus, constructed from data from the Viking and 
Pioneer Venus orbiter mass spectrometers (Nier and McElroy, 1976; Hedin et 
u/., 1983). Above the homopauses, which are located near 120-125 km and 
135-140 km on the Martian and Venusian daysides, respectively, the lighter 
constituents become more important, and atomic oxygen eventually becomes 
the dominant constituent. The Viking mass spectrometers show a transition 



Figure 1.17 Altitude profiles of the major neutral species in the Martian atmosphere. 
The points are from the Viking 1 mass spectrometer data (Nier and McElroy. 1976). 
(Adapted from Fox and Dalgarno, 1979) 
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(a) Altitude (km) 



(b) Altitude (km) 


Figure 1.18 Altitude profiles of the major neutral species in the Venus thermosphere 
for the equator at (a) noon and (b) midnight from a model based on Pioneer Venus 
neutral mass spectrometer data. (Adapted from Hedin et al . , 1983) 


altitude of about 195 km on Mars (Nier and McElroy, 1976, 1977); models 
derived from Pioneer Venus data place the transition altitude near 155 and 
140 km on the dayside and nightside of Venus, respectively (Hedin et aL 1983). 

Early investigations of the ionosphere of Mars predicted that the ion 
composition would be simpler than that of the Earth. The model of McElroy 
(1973b) predicted that, below 250 km, only 0 2 + and C0 2 + would be important; 
the computed densities of N + , NO + , CO + and 0 + were smaller at their maxima 
by one to three orders of magnitude. Only the major ion densities could be 
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derived from data recorded by the retarding potential analyzers on the Viking 
entry probes. Figure 1.19 shows the profiles derived from Viking 1 data (Hanson 
et al . , 1977). Among other instruments, the Pioneer Venus orbiter carried an 
ion mass spectrometer to Venus in 1978. Data were taken inside the ionosphere 
for more than 600 orbits at the beginning of the mission, and again for several 
orbits in the re-entry phase in 1992. An example of the ion density profiles 
obtained on the dayside is shown in Figure 1.20 (Keating et al ., 1985). The ion 
composition is complex. The increased knowledge of ion chemistry on Venus 
obtained from the Pioneer Venus measurements and subsequent modeling 
allowed the ion densities on Mars to be modeled with greater confidence (Nagy 
et al 1980; Fox, 1982a, 1982b, 1985, 1989, 1992b). Except where noted, the 
following discussion refers to the low solar activity model shown in Figure 1.21 
for Mars and the high solar activity model shown in Figure 1.22 for Venus. 



Figure 1.19 Altitude profiles of the ion densities derived from the Viking 1 retarding 
potential analyzer. The filled circles are 0 2 + , the open triangles are C0 2 + and the open 
squares are O*. The heavy solid curve labeled N e is the sum of the individual ion 
densities. The dashed lines through the 0 2 + and C0 2 + data are linear fits to the topside 
data, and exhibit scale heights of 29 and 23 km, respectively. (Adapted from Hanson et 
al . , 1977) 
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Figure 1.20 Ion density profiles for species measured with the Pioneer Venus orbiter 
ion mass spectrometer. The original data have been smoothed. (Adapted from Bauer 
etai, 1985) 



Figure 1.21 Computed altitude profiles of the ions in the ionosphere of Mars for a low 
solar activity model based on Viking 1 data. The neutral densities employed are those 
shown in Figure 1.17 and the solar zenith angle is 45 . The dashed curves are the ion 
densities inferred from the Viking retarding potential analyzer data. (Hanson et al 
1977) 
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Figure 1.22 Computed altitude profiles of ions in the Venus ionosphere for high solar 
activity and 45°. The neutral model employed was that of Hedin el al. (1983) for 
F l01 = 200 


These models were based on the model thermospheres shown in Figures 1.17 
and 1.18. 

On both planets, C0 2 + is produced mostly by photoionization and photo- 
electron impact ionization of C0 2 : 

C0 2 + (ftv, e) ♦ CCV +(e,2e) (153) 

with a peak production rate near 130 km for Mars and 140 km for Venus. On 
Venus, the reaction 

CO + +C0 2 >OV +CO (154) 

exceeds direct production of C0 2 above about 165 km, but on Mars, where 
the CO mixing ratio is smaller, direct production of C0 2 + exceeds production 
by all chemical reactions up to about 205 km, and reaction (154) provides only 
about 10% of the total source at 200 km. It is less important there than charge 
transfer from N 2 * 

N 2 + + C0 2 ► C0 2 + + N 2 (155) 

Reaction (155) is of minor importance in the Venus ionosphere, providing only 
6% of the source from the ion peak to 175 km. At higher altitudes, charge 
transfer from 0 + ( 2 D) and CC( 2 P) 

0 + ( 2 D, 2 P) + C0 2 f O 4- C0 2 + (156) 

is the most important source of C0 2 + . Reaction of 0 + ( 2 D) and 0*( 2 P) with 
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C0 2 can also proceed by 

0 + ( 2 D, 2 P) + C0 2 » 0 2 + + CO (157) 

The rate coefficient for reactions (156) and (157) has been measured (Viggiano 
et al , 1990) to have a value of 1.06 x 10 9 cm 3 s \ with the reaction proceed- 
ing via charge transfer (channel 156) 94% of the time. On Mars, the smaller 
densities of 0 + resulting from the smaller mixing ratio of O renders reactions 
(156) less important; at 200 km, they are comparable to reaction (155). 

In the presence of a small amount of O, C0 2 + is transformed by the reaction 

C0 2 + TO >0 2 + + CO (158) 

but charge transfer may also take place (Fehsenfeld et al , 1970) 

C0 2 + + O >C0 2 + 0 + (159) 

Reaction with O is the most important loss process for C0 2 + over most of the 
photochemical equilibrium regions, which for C0 2 + is below approximately 
235 km on Venus and 210 km on Mars. At higher altitudes, in addition to 

chemical production and loss, transport by diffusion becomes important in 

determining the density profiles. Dissociative recombination 

C0 2 + 4-e >CO + 0 (160) 

is the most effective chemical loss process for C0 2 + above about 225 km on 
Venus. On Mars, although the O mixing ratio is smaller, the total electron 
densities are less by more than an order of magnitude over most of the topside 
ionosphere, and dissociative recombination exceeds reactions with neutral 
species only about 375 km. 

0 + can be produced both by photodissociative ionization of C0 2 and CO 


C0 2 + hv — 

-CO + 0 + +e 

(161) 

CO + hv 

-C + 0 + + e 

(162) 

and by direct photoionization of O 

O + hv — 

— * 0 + +e 

(163) 


On Mars, reaction (159) is more important than direct ionization as a source 
of 0 + below about 180 km. Because of the larger abundance of O in the 
atmosphere of Venus, direct production and reaction (159) are comparable near 
145 km, but at higher altitudes direct production dominates. In addition, the 
reaction 

CO + + O ► 0 + + CO (164) 

also occurring in the supernova mantle, exceeds reaction (159) as a source of 
0 + above 170 km on Venus. Production of 0 + is followed nearly always by 

o 2 + + CO 


0 + + COi 


(165) 
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Reactions (158) and (165) are comparable sources of 0 2 + near the ion peaks 
of both Venus and Mars. At high altitudes, above 190 and 150 km for Mars 
and Venus, respectively, reaction (165) becomes more important. Loss is mainly 
by dissociative recombination 

0 2 + + e ^O + O (166) 

at all altitudes in the photochemical equilibrium region. 

0 2 + may react with both N and NO to produce NO + 

0 2 + + N * NO + 4- O (167) 

0 2 + 4- NO * NO + 4- 0 2 (168) 

Reaction (167) is the largest source of NO + below about 180 km on Mars and 
195 km on Venus. Above those altitudes the reaction 

N 2 + +0 ► NO + 4- N (169) 

is more important. Since NO has the lowest ionization potential of any major 
or minor species in the thermospheres of the terrestrial planets, its loss is almost 
completely by dissociative recombination 

NO + 4e ► N + O (170.) 


N 2 has the very high ionization potential of 15.58 eV, so it can be produced 
chemically in only a few reactions, including 


Ar + + N 2 — 

-N 2 + 

+ Ar 

(171) 

He + + isl- 

-n 2 + 

+ He 

(172) 

and 




0 + ( 2 D) + N 2 — 

-n 2 + 

+ o 

(173) 

Reaction (171) is insignificant as a source of N 2 + 

for Mars, although the 

reverse 


reaction, charge transfer from vibrationally excited N 2 is the major source of 
Ar + at high altitudes. The mixing ratio of Ar on Venus is much smaller than 
that on Mars (see Table 1.2). The densities He and He + have been measured 
by the mass spectrometers on the Pioneer Venus orbiter. It appears that the 
source of N 2 + from reaction (172) is insignificant on the dayside, but it plays 
a larger role in the predawn sector on the nightside of Venus, where the densities 
of light ions maximize. The mixing ratio of He on Mars is unknown and the 
contribution of this reaction is uncertain. 

The rate coefficient for reaction (173) has been measured to be close to the 
gas kinetic value, about 8 x 10 10 cm 3 s _1 (Johnsen and Biondi, 1980; Rowe 
et al , 1980). Reaction (173) is a potentially important source of N 2 + , especially 
at high altitudes where 0 + is produced mainly by photoionization and 
electron-impact ionization of O. Photoionization of O by unattenuated sunlight 
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produces 37% 0 + ( 2 D) and 21% 0 + ( 2 P). Reaction (173) is significant, but not 
dominant at high altitudes on Mars; it is the source of 35% of the N 2 + at 
200 km. On Venus, where the O mixing ratio is higher, charge transfer from 
0 + ( 2 D) is the source of 10% of the N 2 + at the ion peak and the most important 
source above about 175 km. Recently, that conclusion has become more secure 
since the rate coefficients for reactions (156) and (157) have been measured. 
They are the dominant loss processes for 0 + ( 2 D) over most of the ionosphere 
of Mars and below 165 km on Venus. Above that altitude reaction with CO 

0 + ( 2 D) + CO ► CO + + O (174) 

becomes more important. At high altitudes on both planets, radiation and 
quenching of 0 + ( 2 P) are significant sources of 0 + ( 2 D). Some uncertainty still 
exists in the 0 + ( 2 D) budget because the rate coefficients for O + ( 2 D) and O + ( 2 P) 
in reactions (156) and (157) have not been measured separately and because 
the rate coefficient for quenching of 0 + ( 2 D) by O 

O + ( 2 D) + O > 0 + ( 4 S) + O (175) 

has not been measured. A value of 1 x 10 11 cm 3 s' 1 was derived from one 
analysis of Atmosphere Explorer data (Orsini et al , 1977), but another (Oppen- 
heimer et al. , 1976) suggested a larger value of 1 x 10“ 10 cm 3 s [ . This reaction 
is of little importance for the Martian ionosphere, where the densities of O are 
small, but is significant on Venus. At 205 km, reaction (175) comprises either 8 
or 45% the 0 + ( 2 D) loss, for the smaller and larger rate coefficient, respectively. 

The analog of reaction (173) for 0 + ( 2 P) 

0 + ( 2 P) + N 2 - * N 2 + + O (176) 

is also potentially important, but the rate coefficient has not been measured at 
thermal energies. The measured rate coefficient for a mixture of 0 + ( 2 D) and 
0 + ( 2 P) with N 2 is 1.35 x 10“ lo cm 3 s“ l (Glosik et al., 1978). Cross sections 
for the reactions of 0 + ( 4 S), 0 + ( 2 D) and 0 + ( 2 P) with N 2 have been measured 
for relative energies from 5 to 20 eV (Lavollee and Henri, 1989), with the result 
that the ratio of the cross sections for reactions (176) and (173) is about 0.6. 
The importance of reaction (176) as a source of N 2 + , however, cannot be 
accurately ascertained until the rate coefficients for the reactions of 0 + ( 2 P) 
alone with C0 2 and with O have been determined. Measurements of the rate 
coefficients for reactions (175) and (176) would greatly improve our under- 
standing of the N 2 + densities in the ionospheres of the Earth and Venus. 

N 2 + is destroyed mainly by charge transfer to C0 2 (reaction 155) in the 
Martian ionosphere, but reaction with O (reaction 169) is also significant at 
high altitudes, providing more than 20% of the total sink above 200 km. The 
chemistry of N 2 + is more complicated on Venus. At low altitudes reaction ( 1 55) 
is the dominant sink for N 2 4 , but above about 170 km, reaction (169) becomes 
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more important. Reaction with CO 

N 2 + +CO >CO + +N 2 (177) 

is a minor loss process for N 2 + on Venus. Dissociative recombination 

N 2 + + e ► N + N (178) 

is insignificant as an N 2 + sink in the ionosphere of Mars, for the same reasons 
previously given for dissociative recombination of C0 2 + , but is comparable to 
reaction (169) above 200 km on Venus. 

The ionization potential of CO is about 14eV. On Venus, the production of 
CO + is dominated by direct photoionization and electron-impact ionization of 
CO, but the reaction 

C + + C0 2 * CO + + CO (179) 

also contributes significantly, making up about 25% of the total source at the 
ion peak and somewhat less at higher altitudes. The major source above about 
175 km is charge transfer from 0 + ( 2 D) to CO: 

0 + ( 2 D) + CO ► CCU + O (180) 

On Mars, where the mixing ratios of O and CO are smaller, reaction (179) 
provides a larger fraction, about 40% of the production of CO + at the ion 
peak, and reaction (180) is negligible as a source of CO + . 

The contribution of the reaction of atomic carbon with 0 2 + 

C + 0 2 + >CO + + O (181) 

is a potentially important source of CO + in the ionosphere of Venus, especially 
in the region of the ion peak, but the rate coefficient is unknown and the possible 
products include C + as well as CO + . Using the estimated rate coefficient of 
5 x 10 -n cm 3 s“ 1 (Prasad and Huntress, 1980), we find that reaction (181) 
produces 5% of the CO + at 145 km. If the reaction actually proceeds at near 
gas kinetic rates, about 1 x 10~ 9 cm 3 s“ \ the source would be much larger. 
Loss of CO + is by reaction with C0 2 

CO + + C0 2 >C0 2 + + CO (182) 

at the ion peaks of both Mars and Venus, but charge transfer to O becomes 
more important above about 175 km on Venus. On Mars, reaction (164) is a 
minor sink, providing more than 20% of the total loss only above 200 km. 
Dissociative recombination (Mitchell, 1990) 

CO + + e >C + 0 (183) 

is not a significant loss process at any altitude in the Martian ionosphere, but 
is the dominant loss process for CO + above about 250 km on Venus. 

Density profiles of C + were measured in situ by the ion mass spectrometer 
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on the Pioneer Venus orbiter (Taylor et al , 1980), but no measurements have 
been carried out in the Mars ionosphere. C + is produced in dissociative 
photoionization of C0 2 and CO and in direct photoionization and electron- 
impact ionization of atomic carbon. The atomic carbon densities on Venus and 
Mars have not been measured in situ , but the ultraviolet emission lines of C 
were observed by the Mariner and Pioneer Venus spacecraft (Barth et al ., 1971 ; 
Barth et al . , 1972; Paxton, 1985) and in ultraviolet spectra of Venus taken from 
a rocket (Rottman and Moos, 1973). The major source of C in the Venus 
thermosphere and at high altitudes in the Martian thermosphere is photo- 
dissociation of CO 


CO + hv - 

-^c + o 

(184) 

but photodissociation of C0 2 



CO 2 -I- hv 

- c + o + o 

(185) 


also contributes, especially near the ion peak on Mars (Krasnopolsky, 1982, 
1983; Fox, 1982b). 

The major source of C + is direct ionization at all altitudes on Mars. On 
Venus the reaction 

O/ +C ^0 2 + C + (186) 

provides 25% of the total source near the ion peak and is comparable to direct 
production near 175 km if the rate coefficient is near the estimated value of 
5 x l(T n cm 3 s -1 (Prasad and Huntress, 1980). Photoionization of atomic 
carbon is the most important production mechanism for C + above about 
200 km (Fox, 1982b). There is also a potential source due to the reaction 
(Paxton, 1983) 

0 + +C >C + -bO (187) 

The process may also be an important source of C + in the supernova, but the 
rate coefficient is unknown. The reaction will dominate the production of C + 
at high altitudes on Venus if the rate coefficient is greater than 1 x 10“ 10 
cm 3 s _1 , but will be insignificant if the rate coefficient is less than 1 x 10“ 11 
cm 3 s“ V The C + budgets on Venus and Mars cannot be completely understood 
until the rate coefficients for reactions (186) and (187) have been measured. 

Density profiles of N + with peak values sometimes exceeding 10 3 cm -3 have 
been measured in the Venus ionosphere (Taylor et al ., 1980). The major sources 
below about 200 km are dissociative ionization of N 2 

N 2 + (/iv,e) ► N + + N +(e,2e) (188) 

and photoionization and electron-impact ionization of N 
N + (/? v, e) ► N + -f (e, 2e) 


( 189 ) 
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Chemical production of N + in the reactions 


0 2 + + N( 2 D) — 

— » N + + 0 2 

(190) 

0 + ( 2 D)+ N — 

— N + + O 

(191) 

0 + ( 2 P) + n 2 — 

— ► N + + NO 

(192) 


provides about 14% of the source at the ion peak and exceeds direct production 
above about 200 km. Reaction (190) is nearly resonant and the rate coefficient 
can be derived from the measured value of the reverse reaction (Dalgarno, 
1970); a value of 9.3 x 10 11 cm 3 s" 1 is implied by the measurements (Smith 
et al . , 1983b; O’Keefe et al . , 1986). The rate coefficient for reactions (191) and 
(192) have not been measured, and values estimated from Atmosphere Explorer 
Data (Dalgarno, 1979) are used in the model. Measurements of these rate 
coefficients would be valuable to interpreting the measured N + densities for 
both Venus and Earth, and predicting the densities for Mars. 

Loss of N + on Venus is dominated at low altitudes by reaction with C0 2 , 
which proceeds at nearly gas kinetic rates by 


N + + CC^- 

— ♦ N + C0 2 + 

(193) 

and 



N + + C0 2 — 

— NO + +CO 

(194) 

The branching ratios between reactions (193) and (194) are 75 and 
respectively (Smith et al . , 1978). At altitudes above 175 km, reaction of N 
CO 

25%, 
+ with 

N + + CO — 

— CO + + N 

(195a) 

— 

-* NO + + C 

(195b) 


is more important than that with C0 2 . Reaction (195b) is also important in 
dense interstellar cloud chemistry. Reaction with O 

N + 4- O * N + O 1 (196) 

is a minor loss mechanism, providing only 5% of the total N f sink at 200 km. 

The chemistry of N + in the Martian ionosphere is much simpler due in part 
to the lower thermospheric mixing ratios of O and CO. Its production is largely 
by ionization of N (reaction 188) and dissociative ionization of N 2 (reaction 
189). The only significant chemical source is reaction of 0 2 + with N( 2 D) 
(reaction 190). Loss is mainly by reaction with C0 2 (reactions 193 and 194). 
Reactions (195a) and ( 1 95b) are not important, accounting for only 5% of the 
loss at 200 km. 
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1.4.2. 2 Nightside ionospheres 

The nightside ionosphere of Venus was first detected by the radio occultation 
experiments on the Mariners 5 and 10 spacecraft (Kliore et al , 1967; Fjeldbo 
et al ., 1975) and more recently by the Venera 9 and 10 spacecraft and the 
Pioneer Venus orbiter (Gringauz et al ., 1977; Kliore et al ., 1979). In situ 
measurements of the ion densities were carried out by the PV orbiter ion mass 
spectrometer (OIMS) (Taylor et al , 1980). Although the same constituents were 
detected by the OIMS on the nightside as on the dayside, the densities of heavy 
atomic ions, such as O*, N + and C + , were observed to be about an order of 
magnitude less, and the molecular ion densities about two orders of magnitude 
less than dayside values. 

The densities of light atomic ions, such as He* and H + , were observed to 
maximize in the predawn sector. Both the neutral densities, as measured by the 
neutral mass spectrometer, and the ion densities on the nightside were observed 
to be highly variable (Niemann et al , 1980). 

The nightside ionosphere of Venus is produced by some combination of 
plasma transport from the dayside (Knudsen et al. , 1980; Spenner et al , 1981; 
Cravens et al , 1983) and precipitation of electrons that have been observed in 
the umbra (Chen and Nagy, 1978; Gringauz et al , 1979; Breus et al , 1985; 
Knudsen and Miller, 1985). Probably the former is more important at high 
solar activity and the latter at a low solar activity (Knudsen, 1988; Kliore et 
al , 1991). Even at high solar activity, however, the source due to electron 
precipitation varies (Fox and Taylor, 1990). Fox (1992b) has discussed the 
chemistry of the nightside ionosphere that is produced by transport of atomic 
ions from the dayside. In this case, the ion sources are chemical. Ions that are 
produced from parents with high ionization potentials, such as N 2 + , CO* and 
C0 2 + , are difficult to produce chemically. For example, in the absence of 
ultraviolet photons or energetic electrons, the only potential sources of N 2 * are 

He* + N 2 ► N 2 * + He (197) 

O 2 * + N 2 * N 2 * + O* (198) 

and 

0*( 2 D) + N 2 ► N 2 * + O (199) 

The products of reaction ( 1 98) are unknown, but the large exothermicity implies 
that dissociative channels probably dominate. It is also likely that 0*( 2 D) does 
not survive transport to the nightside; consequently reaction (199) probably is 
not important. Only reaction (197) is found to be important, especially in the 
predawn region where the densities of light ions are seen to maximize. 

For CO* there are more chemical production mechanisms, including reac- 



52 


A. DALGARNO AND J. L. FOX 


tions, that are analogous to reactions (197) to (199) 

He f 4- CO z ♦ He + CO + + O 

(200) 

o 2 + 4- CO — 

— co + +o 

(201) 

0 + ( 2 D + ) + CO — 

— CO + 4-0 

(202) 


The reaction of He + with CO proceeds mostly by dissociative charge transfer 
(67) 

He + +CO >C + + O + He 

and is therefore not a source of CO + . The same comments apply to reactions 
(201) and (202) as to reactions (198) and (199). The most important source of 
CO + near the peak of the nightside ionosphere is the reaction 

C 4 + C0 2 * CO 4 + CO (203) 

At high altitudes, the most important source is charge transfer from N + to CO 
in reaction (195a). The fractions of C f and transported to the nightside 
ionosphere are, however, small compared to 0 + . The dearth of chemical 
production mechanisms for CO + and N 2 + and their large variability in the 
nightside OIMS data lead to the suggestion (Fox and Taylor, 1990) that the 
appearance of large densities of mass-28 ions is a signature of auroral precipita- 
tion. The sources of CO 4 and N 2 + for the dayside and nightside ionospheres 
are compared in Figures 1.23(a) to (d) and 1.24(a) to (c). 

The major chemical production mechanisms of C0 2 + on the nightside of 
Venus are charge transfer from N + 

N 4 + C0 2 * C0 2 4 + N (204) 

and from the mass-28 ions CO 4 and N 2 4 

CO 4 + C0 2 ► CO + C0 2 4 (205) 

and 

N 2 + + C0 2 * C0 2 4 4- N 2 (206) 

Reaction (205) dominates near the ion peak, but reaction (204) is more 
important at high altitudes. C0 2 + has the potential to be enhanced substanti- 
ally in the presence of electron precipitation. In the plasma-transport dominated 
ionosphere, 0 2 4 is produced mainly by 

O 4 + C0 2 >0 2 4 +CO (207) 

In the presence of electron precipitation, the reaction (158) is the major source 
in the lower ionosphere. The sources of C0 2 + in the dayside and nightside 
ionospheres are compared in Figures 1.23(e) and (1) and 1.24(d) and (e). 

NO 4 is produced mainly by the reaction 

0 2 4 + N ► NO 4 + O 


( 208 ) 






{Ui>|} apnjmv 


(d) log(N 2 + loss rate){cm' 3 s' 1 } (©) log (C0 2 + production rate) (cm' 3 s' 1 ) (f) log(C0 2 + loss rate ) (cm 3 s' 1 ) 

Figure 1.23 Altitude profiles for production and loss of ions for a high solar activity model of the Venus dayside ionosphere 
(a) production of CO + , (b) loss of CO + , (c) production of N 2 + , (d) loss of N 2 + , (e) production of C0 2 + , (0 loss of C0 2 + 
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under low precipitation conditions. Its loss is always by dissociative recombina- 
tion. Hence its density is given approximately by 

[NO + ] = fc[N][0 2 4 ]/a[e] (209) 

where k is the rate coefficient for reaction (208) and a is the dissociative 
recombination coefficient. Near the ion peak where [e] — [Oj 1 ], the NO + 
density is largely independent of 0 + flux. Computed ion density profiles 
produced by imposing a downward flux of 0 + of 1 x 10 8 crrT 2 s _1 with 
smaller fluxes of C + , N + and He + near the antisolar point are shown in Figure 
1.25. 

The nightside ionosphere of Mars has been detected at low solar activity by 
radio occultation experiments on the Viking and Soviet Mars 4 space probes 
(Vasil’ev et al , 1975; Zhang et ai % 1990). Because the orbit of Mars lies outside 
that of the Earth, only a small solar zenith angle range on the nightside past 
the terminators has been investigated. A peak did not always appear, but when 
it did the peak density averaged about 5 x 10 3 cm -3 . The HARP instrument 
on the Soviet Phobos spacecraft measured large fluxes of electrons in the vicinity 
of Mars (Verigin et al . , 1991). A two-stream code was used to model the 
ionization rates and perform a photochemical equilibrium calculation of the 



Figure 1.25 Computed altitude profiles of the major ions in the Venus nightside 
ionosphere resulting from the transport of atomic ions from the dayside. Downward 
fluxes ofO*,C*, N* and He* were imposed at the upper boundary of the model, near 
235 km. The neutral model was taken from Hedin et al. (1983) and is appropriate to 
the equator at 2300 h local time 
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ion and electron density profiles that would result if the observed electrons 
precipitated into the nightside atmosphere (Haider et ai , 1992). 

Models of the dayside ionosphere of Mars have shown that the ion density 
profiles cannot be reproduced without imposing a loss process at the upper 
boundary of the model (Fox, 1993; Shinagawa and Cravens, 1992). It has been 
proposed that this loss process is the divergence of the horizontal ion fluxes, 
by analogy to Venus (Shinagawa and Cravens, 1989, 1992). In order to compute 
an upper limit to the nightside ion densities. Fox et ai (1993) imposed, in the 
downward direction at the upper boundary of a model of the nightside 
thermosphere, the upward fluxes O 4 , C + , N + and 0 2 + that were found 
necessary (Fox, 1993) to reproduce the dayside ion density profiles. Fox et ai 
(1993) also performed a calculation of the ion densities that would result from 
electron precipitation using a multistream electron transport code. The iono- 
sphere produced by ion transport alone showed a peak density of about 
1.3 x 10 4 cm 3 at an altitude of 179 km. When combined with electron pre- 
cipitation, the peak ion density increased to 1. 8-1.9 x 10 4 cm -3 . Altitude 
profiles of the ion densities for plasma transport alone and for precipitation of 
the electrons detected by the Phobos HARP instrument in the magnetosphere 
lobes are shown in Figure 1.26(a) and (b). 




(b) Martian nightside ion densities (cm' 3 ) 

Figure 1.26 Computed nightside ion densities of Mars (a) due to transport of ions from 
the dayside and (b) precipitation of the electrons detected in the magnetotail lobes by 
the Phobos HARP instrument (Verigin et ai , 1991). (After Fox et ai , 1993) 
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1.4.3 TITAN AND TRITON 

The atmosphere of Titan is composed mostly of N 2 , with small amounts of 
methane and higher hydrocarbons, as shown by the model thermosphere in 
Figure 1.27. In this model, the N 2 and CH 4 density profiles were taken from 
the Voyager solar occultation measurements (Smith et al . , 1982), and those of 
C 2 H 2 , C 2 H 4 and C 2 H 6 were adopted from a photochemical model (Yung et 
al . , 1984). N 2 is the major species below 1700 km and CH 4 is relatively more 
abundant at higher altitudes. In contrast to the thermosphere of Jupiter, the 
mixing ratios of the most important hydrocarbons, methane, acetylene, and 
ethylene do not decrease above the Titan homopause, which is near 925 km. 
Methane is lighter than N 2 , and is therefore more abundant near the altitude 
of peak ion production (about 1050 km) than in the bulk atmosphere. On 
Jupiter, the peak photoionization of H 2 is near 600 km, almost 200 km above 
the homopause. Since methane and other hydrocarbons are heavier than H 2 , 
they are severely depleted in this altitude region and are not involved in the 
ion chemistry near the peak. On Titan, however, hydrocarbons are important 
in the ion chemistry above, near and below the peak. 

Models of the ionosphere of Titan have been constructed (Ip, 1990a; Fox 
and Yelle, 1991; Keller et al ., 1992). Solar photoionization and photoelectron 
impact ionization are sufficient to account for the electron density profiles 
measured by the Voyager radio science experiment (Keller et al ., 1992). The 



Figure 1.27 Neutral model of the thermosphere of Titan. The densities of N 2 and CH 4 
were taken from the solar occultation profiles of Smith et al. (1982). The density profiles 
of ethane, ethylene and acetylene were taken from the photochemical model of Yung et 
al. (1984) 
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log (ton production rate) (cm ’ 3 s' 1 ) 

Figure 1.28 Computed production rate profiles of the major ions in the ionosphere of 
Titan for the neutral thermosphere shown in Figure 1.27. (From Fox and Yelle, 1991) 

production rates of six ions due to photoionization and photoelectron impact 
ionization (Fox and Yelle, 1991) are shown in Figure 1.28. Although N 2 + is 
the dominant ion produced below 1700 km, because the ionization potential of 
N 2 is very large (15.6 eV), N 2 + is not expected to be the terminal ion. Altitude 

profiles of the major loss processes for N 2 + in the ionosphere of Titan are 

shown in Figure 1.29. N 2 + is lost in reaction with CH 4 , largely by 

N 2 + + CH 4 > CH 3 + + N 2 + H (210) 

but also by 

N 2 + + CH 4 ► CH 2 + +N 2 + H 2 (211) 

The reaction of N 2 + with H 2 is also rapid 

N 2 + 4- H 2 ► N 2 H + + H (212) 

Charge transfer to acetylene 

N 2 + + C 2 H 2 >C 2 H 2 + +N 2 (213) 

and to HCN 

N 2 + + HCN ► HCN + + N 2 (214) 

and dissociative recombination 

N 2 + + e * N 4- N 

are smaller but significant destruction mechanisms. 


( 215 ) 
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Figure 1.29 Computed loss rate for N 2 + in the ionosphere of Titan. Ion molecule 
reactions are much more important than dissociative recombination, for which the 
maximum destruction rate is about two orders of magnitude less than reaction with 
methane 

The N 2 H + produced in reaction (212) above may react with methane to 
produce CH 5 + 

N 2 H + + CH 4 >N 2 + CH 5 + (216) 

Figure 1.28 shows that, at high altitudes, CH 4 + is the most important ion 
produced; it also reacts at near gas kinetic rates with CH 4 to produce CH 5 + 


ch 4 + + ch 4 - 

— » ch 5 + +ch 3 

(217) 

As a result of these reactions, CH 5 + is the most important ion at high altitudes 
in the Titan ionosphere (Fox and Yelle, 1991; Keller et a/., 1992). It is lost by 
reaction with acetylene 

ch 5 + + c 2 h, — 

— ♦ c 2 h 3 + + ch 4 

(218) 

and with ethylene 



CH 5 + + C 2 H 4 - 

— ♦ c 2 h 5 + + ch 4 

(219) 

and by dissociative recombination 



CH 5 + + e- 

— ♦ CH 4 + H 

(220a) 

- 

— * CH , + 2H 

(220b) 


The product channels in reactions (220a) and (220b) have not been determined, 
but a yield of 1.19 H atoms per dissociative recombination of CH 5 + has been 
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measured (Adams et al, 1991) and it has been assumed that reactions (220a) 
and (220b) proceed by channel (220a) 81% of the time and by channel (220b) 
19% of the time (Fox and Yelle, 1991). 

Since both cis- and tram - HCNH have very low ionization potentials of 6.8 
and 7.0 eV, respectively (Nesbitt et aL, 1991), there are many production 
mechanisms for HCNH + . The most important pathways leading to the produc- 
tion of HCNH + are shown in Figure 1.30, and altitude profiles of the most 
important production and loss mechanisms are presented in Figure 1.31. The 
most important sources are the proton transfer reactions 


C 2 H 5 + + HCN - 

— * C 2 H 4 + HCNH + 

(221) 

and 



C 2 H 3 + 4- HCN- 

— » C 2 H 2 + HCNH + 

(222) 

Of somewhat less importance are the reactions 


N + + CH 4 — 

— » HCNH + + CH 2 

(223) 

and 



HCN + + CH 4 — 

— > HCNH + + CH 3 

(224) 

Dissociative recombination is the only significant loss mechanism over the 
entire ionosphere. There are four exothermic product channels for this reaction : 

HCNH + + e 

► HCN + H 

(225a) 


♦ HNC + H 

(225b) 


► CH, + N 

(225c) 


> H 2 + CN (22 5d) 

Measurements (Adams et ai , 1991) indicate that 0.63 H atoms are produced 
per recombination. This implies that the sum of the branching ratios for 
channels (225a) and (225b) is 63% and the sum of the branching ratios for 
channels (225c) and (225d) is 37%. One model (Fox and Yelle, 1991) was 
constructed with the arbitrary assumption that production of HCN, the more 
stable isomer, is the dominant channel of the first two and that the third channel 
is more important than the fourth, whereas another (Keller et al . , 1992) assumed 
that all the reactions proceed via channel (225a). 

CH 3 + is produced by dissociative ionization of methane and by photoioniza- 
tion of CH 3 by the strong solar Lyman alpha line. It reacts rapidly with methane 
to produce C 2 H 5 + 

CH 3 + + CH 4 ► C 2 H 5 + + H 2 (226) 

C 2 H 5 + is also produced in proton transfer reactions to ethylene, the most 
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log (production and loss rates)(crrf 3 s* 1 ) 

Figure 1.31 Altitude profiles for the major production and loss mechanisms for 
HCNH + . The solid curves are production mechanisms and the dashed curve is for the 
only loss process, dissociative recombination 


important of which are 

CH 5 + + C 2 H 4 >C 2 H 5 + + CH 4 (227) 

and 

C 2 H 3 + + C 2 H 4 >C 2 H 5 + +C 2 H 2 . (228) 

C 2 H 5 + is destroyed mainly in reactions with HCN (reaction 222) with acetylene 

C 2 H 5 + -F C 2 H 2 * C 2 H 3 + + C 2 H 4 (229) 

and in dissociative recombination 

C 2 H 5 + + e > C 2 H 4 + H (230a) 

► CH 3 + CH 2 (230b) 


Here, again, the actual distribution of the products in dissociative recombina- 
tion is uncertain. In one model (Fox and Yelle, 1991), the two channels shown 
above are assumed to be of equal importance and in another (Keller et al , 
1992), the reactions are assumed to produce C 2 H 2 + H 2 (or 2H) + H. 

Altitude profiles of the eight most important ions predicted for the Titan 
ionosphere at high solar activity are shown in Figure 1.32 (Fox and Yelle, 1991). 
The most abundant ions are H 2 CN + below about 1200 km and CH 5 + above 
about 1200 km. The peak electron density is about 6 x 10 3 cm -3 . Upper limits 
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Figure 1.32 Computed density profiles of the major ions in the ionosphere of Titan. 
The solid curves are ion densities and the dashed curve is the total electron density 
profile. (From Fox and Yelle, 1991) 


to the electron density at the morning and evening terminators of 5 x 10 3 and 
3 x 10 3 cm -3 , respectively, were obtained by the Voyager 1 radio occultation 
experiment (Lindal et al . , 1983). The model shown in Figure 1.32 is for a solar 
zenith angle of 60 and is consistent with the Voyager upper limits. 

The atmosphere of Triton is similar to that of Titan in that the major 
constituent is N 2 . In N 2 atmospheres, usually the minor constituents control 
the absorption of radiation. For Earth the absorbing species is 0 2 and for 
Titan, CH 4 . Triton, however, unlike Titan or Earth, appears to have an extreme 
dearth of minor constituents. A model of the Triton atmosphere is shown in 
Figure 1.33 (Krasnopolsky et al 1993). The surface pressure is only about 
14 pbar (Gurrola et al . , 1991), and methane in the troposphere is in equilibrium 
with a surface methane frost at about 38 K. The tropopause temperature is 
only about 2 K colder than the surface and occurs in the altitude range 8-12 km. 
The base of the thermosphere is near the tropopause, and temperatures increase 
monotonically above the tropopause to an exospheric temperature of about 
100 K. The methane profiles derived from Voyager stellar occultation data 
(Herbert and Sandel, 1991) indicate that the methane density decreases more 
rapidly than the N 2 density from a surface value near 5 x 10 1 1 cm~ 3 to about 
10 6 cm 3 at 80 km. The higher hydrocarbons, acetylene, ethylene and ethane, 
condense out between 60 and 90 km. H 2 and H are produced by Lyman alpha 
photolysis of methane; they diffuse upward toward the exobase near 800- 
900 km, where they may escape thermally. 
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Figure 1.33 Model of the Triton atmosphere based on Voyager Ultraviolet Spectro- 
meter measurements. The CO mixing ratio is assumed to be 10" 3 . The N 2 , N and 
temperatures profiles are from solar occultation measurements, and the H and H 2 
densities are from calculations. The short dashed curves are for solar ionization only, 
while the solid and long-dashed curves are for different assumptions about the interac- 
tion of magnetospheric electrons with Triton’s thermosphere. The solid curves are the 
recommended profiles. (Adapted from Krasnopolsky et al . . 1993) 


On Triton, the principal ions produced are N 2 + and N + . On Titan, N + is 
lost primarily by the reaction with CH 4 , via reaction (223) or by 

N" + CH 4 ► CH 3 + + NH (231a) 

>HCH + +H 2 + H (231b) 

► CH 4 + + N (231c) 

Most of the reactions proceed by channel (223) or (231a) (Anicich et al , 1977). 
The reaction with C 2 H 4 is a minor loss process, as is the reaction 

N + +H 2 ► NH + + H (232) 

Radiative recombination 

N + + e >N + /iv (233) 

is very slow. On Triton, the CH 4 mixing ratio is too small for reactions (223) 
and (231) to contribute significantly to the loss of N + , and reaction (232) with 
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H 2 is more important. N 2 + is also destroyed by reaction with H 2 : 
N 2 + + H 2 * N 2 H + + H 


with H 


(234) 


N 2 + + H ► N 2 4- H + (235) 

and in dissociative recombination (215). Because the N + loss rates are very 
small, N + has been predicted to be the major ion near the ionospheric peak 
of Triton in several models (Ip, 1990b; Majeed et al ., 1990; Yung and Lyons, 
1990). 

It appears that, unlike Titan, solar photoionization and photoelectron impact 
ionization are not sufficient to produce the electron density profiles inferred 
from the Voyager radio science experiments, and ionization by energetic 
electrons from Neptune's magnetosphere has been suggested as the missing 
source (Tyler et at., 1989). The ingress and egress profiles exhibit maximum 
electron densities near 350 km of 2.3 x 10 4 and 4.6 x 10 4 cm -3 , respectively. 
The model shown in Figure 1.34 was constructed by varying the ionization rate 
to fit the Voyager radio science ingress profile; a column ionization rate (in 
addition to the solar source) of 3 x 10 8 (Majeed et al ., 1990) or 4 x 10 8 
cm” 2 s" 1 was inferred (Strobel and Summers, 1993). In order to fit the topside 
scale height of the electron density profile, it was necessary to impose a loss 



Figure 1.34 Model ionosphere of Triton constructed assuming an ionization rate due 
to magnetospheric electrons of 3 x 10 cm 2 s _l , in addition to the solar source. The 
curve labeled RSS is the electron density profile obtained from the Voyager radio 
occultation experiment (Herbert and Sandel, 1991). (From Majeed et al 1990) 
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process at the top of the ionosphere of 4.3 x 10 7 N + ions cm 2 s 1 (Yung and 
Lyons, 1990). Cluster ions of the form N 2 + (N 2 )„ and N + (N 2 )„, which are 
strongly favored energetically, should readily form at the low temperatures of 
the Triton ionosphere. Such clusters may make up most of the ions below about 
200 km (Delitsky et ai , 1990). However, the chemistry of trace species such as 
H 2 , CH 4 and CO modify the conclusion, which may be valid only in the 
troposphere (Strobel and Summers, 1993). 

The Triton ionosphere could be maintained by solar radiation alone if the 
major ion was C + , produced by photoionization of C and by charge transfer 
from N 2 + (Lyons et al . , 1990) 

N 2 + + C ► C + + N 2 (236) 

The rate coefficient for reaction (236) has not been measured; a near gas kinetic 
value is likely. The atomic carbon in the model may originate ultimately from 
photolysis of methane below 100 km, and be destroyed by photoionization, by 
ion-molecule reactions, such as reaction (236), and by upward diffusion and 
Jeans (thermal) escape near the exobase. The major loss process for C + is 
reaction with methane 


c + +ch 4 

- C 2 H 3 + + H 

(237a) 

— 

- c 2 h 2 + + h 2 

(237b) 

The computed C + profile exhibited 
250 km. The addition of the reaction 

a peak density of 3 x 10 4 cm 

3 near 

C + N 2 4- N 2 - 

♦ CNN + N 2 

(238) 


to the model reduced the computed atomic carbon densities, and thus the C + 
densities also, by an order of magnitude (Strobel and Summers, 1993). 


1.4.4 THE EARTH 

The division of the terrestrial atmosphere into regions is based on the structure 
in the altitude profile of the electron density, which is produced by partially 
overlapping layers of ions. These layers arise from altitude variations both in 
the composition of the thermosphere and in the relative importance of different 
ionization sources. Figure 1.35 shows a schematic diagram of the major 
terrestrial ionospheric regions and their most important ionization sources 
(Bauer, 1986). In the lower layers, below about 300 km, the ionosphere is 
essentially in photochemical equilibrium, and transport processes do not affect 
substantially the ion and electron density profiles. Molecular ions are the 
dominant constituents of the D, E and F, regions. The rate of plasma transport 
by ambipolar diffusion increases with decreasing density and, at sufficiently 
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Figure 1.35 Schematic diagram of the regions of the terrestrial ionosphere and their 
principal ionization sources. (After Bauer, 1973) 

high altitudes, the major loss process for atomic ions is diffusion. This results 
in the formation of a peak in the topmost layer of the ionosphere, the F 2 region. 
Here the major ion is 0 + , and the peak electron density can exceed 10 6 cm -3 . 
A model of the F, and F 2 regions of the ionosphere at low solar activity from 
100 to 450 km is shown in Figure 1.36. Most of the discussion of sources and 



log (ion Densities)(cm 3 ) 

Figure 1.36 Altitude profiles of the densities of the ions in the terrestrial F, and F 2 
regions for a low solar activity model of the terrestrial ionosphere 
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sinks that follows is based on this model. At higher solar activity, the neutral 
temperatures are larger, resulting in a more extended thermosphere and 
ionosphere than that shown in Figure 1.36. Also, data from ground-based 
experiments and spacecraft show considerable spatial and temporal variability 
in the ion and electron density profiles. 

NO + is the most abundant ion in the region of the ionosphere. Altitude 
profiles of its sources and sinks are shown in Figure 1.37. The reactions 


n 2 + + o- 

► NO + + N 

(239) 

o + + n 2 - 

♦ NO + + N 

(240) 


produce more than 90% of the NO + above 200 km (Breig et al , 1984). In the 
low solar activity model shown in Figure 1.36, reaction (239) is the most 
important production mechanism below about 275 km, and at higher altitudes, 
reaction (240) dominates. Below about 150 km, a number of reactions become 
significant, including charge transfer from 0 2 + to NO (reaction 168), the 
reaction 

N * 4- 0 2 ► NO + 4 O (241) 

and the reaction of 0 2 + with N (reaction 208), which is the major source of 
NO + at low altitudes on Mars and Venus. Near 100 km, the reaction 

0 2 + 4 N 2 * NO + 4 NO (242) 

may also become important if the rate coefficient is near the upper limit of 
1 x 10~ 15 cm 3 s -1 reported by Ferguson et al. (1965). As in the ionospheres 



Figure 1.37 Altitude profiles of the production (solid curves) and loss (dashed curve) 
of NO + or a low solar activity model of the terrestrial dayside ionosphere 
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of Venus and Mars, loss of NO + is predominantly by dissociative recombina- 
tion at all altitudes. 

The primary sources and sinks of 0 2 + for a low solar activity model are 
shown in Figure 1.38. 0 2 + is produced mainly by photoionization and photo- 
electron impact ionization of 0 2 below about 150 km. At higher altitudes, the 
reaction 

O* + 0 2 * 0 2 + + O (243) 

is more important. Reaction (243) is exothermic by 1.55 eV, and 0 2 + vibrational 
levels up to v = 7 are energetically accessible. 0 2 + has been shown to be 
significantly vibrationally excited at altitudes above about 250 km (Fox, 1986a). 
Reaction (243) is also the major source of 0 2 + (v = 1) in the altitude range of 
about 170-325 km. The reaction 

N + + 0 2 ► N( 2 D, 4 S) + 0 2 + (244) 

can produce vibrationally excited 0 2 + only when the N atom is produced in 
the 4 S ground state. Reaction (244) is the most important source of 0 2 + (r = 1) 




log 0 2 + Loss Rate (cm' 3 s' 1 ) 


Figure 1.38 Altitude profiles of the (a) sources and (b) sinks of 0 2 + for a low solar 
activity model of the terrestrial dayside ionosphere 
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above 325 km and of total 0 2 + above about 450 km. The reaction 

0 + ( 2 D) 4- 0 2 > 0 2 + + O (245) 

has been found to proceed at near gas kinetic rates (Rowe et a /., 1980), but is 
only a minor source of 0 2 + at high altitudes, producing 5% of the 0 2 + at 
250 km and about 8% near 350 km. Charge transfer from N 2 + to 0 2 

N 2 + +0 2 >0 2 + +N 2 (246) 

is significant only at low altitudes. Surprisingly, the reaction of 0 + with C0 2 
(reaction 165) is nearly as important as reaction (246) near 150 km, where it 
constitutes about 8% of the total source of 0 2 + . 

Loss of 0 2 + is mainly by dissociative recombination (reaction 166) at all 
altitudes, but the reactions with NO and N (reactions 167 and 168) together 
contribute about 15% of the total sink near 100 km. Reaction with N 2 (reaction 
242) may also be important at low altitudes. 

Altitude profiles of the primary sources and sinks of N 2 + are shown in Figure 
1.39. Below about 250 km, most of the N 2 + is produced by direct photoioniza- 




log N 2 + Loss Rate (cm * 3 s' 1 ) 

Figure 1.39 Altitude profiles of the primary (a) sources and (b) sinks of N 2 + for a low 
solar activity model of the terrestrial dayside ionosphere 
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tion and photoelectron impact ionization of N 2 . Charge transfer from 0 + ( 2 D) 

0 + ( 2 D) + N 2 >N 2 + +0 (247) 

becomes more important at higher altitudes. If the rate coefficient measured 
recently in the 8-20 eV range (Lavollee and Henri, 1989) also applies to thermal 
energies, the analogous charge transfer reaction from 0 + ( 2 P) 

0 + ( 2 P) + N 2 * N 2 + + O (248) 

provides about 5% of the total production between 150 and 250 km. The 
sources due to charge transfer from He 1 (reaction 172) and from Ar + (reaction 
171) are insignificant. The reverse of reaction (171), however, is exothermic for 
N 2 + in vibrational levels r > 0, and has been found to be the most important 
source of Ar + at high altitudes (Fox, 1986a). 

Loss of N 2 4 is mostly by reaction with O 

N 2 + + O * NO + + N (249a) 

► N 2 + CC (249b) 

up to about 280 km. The major product channel is (249a); only about 10% of 
the reactions proceed via charge transfer (Knutsen et a/., 1988). At higher 
altitudes dissociative recombination (reaction 1 78) is the dominant loss process. 
Near 100 km, charge transfer to 0 2 (reaction 246) becomes significant, provi- 
ding about 10% of the total sink of N 2 + . 

Analyses based on Atmospheric Explorer (AE) data showed that model N 2 + 
densities were larger than the measured values by a factor of about 2 (Breig et 
ai , 1983; Abdou et ai , 1984). Suggestions for improving the agreement included 
increasing the dissociative recombination coefficient for N 2 + (about 2 x 10” 7 
(300/7^) 0,39 cm 3 s _l ) (Zipf, 1980) by a factor of 2-3 and enhancing the rate 
coefficient for charge transfer to O (reaction 249b). Either end could be at least 
partially achieved with accepted values for the rate coefficients of ground state 
N 2 + if vibrationally excited ions reacted much faster than ground state ions. 
Fox and Dalgarno (1985) modeled the vibrational distribution of N 2 + and 
showed, however, that enhanced loss of vibrationally excited N 2 + could reduce 
the densities of N 2 + by a maximum of about 30%. Recently a slightly larger 
dissociative recombination coefficient for ground state N 2 + of about 2.6 x 10” 7 
cm 3 s _1 at 300 K was measured (Canosa et ai , 1991) but a large ab initio 
calculation gives a smaller value of about 1.6 x 10” 7 cm 3 s” 1 (Guberman, 
1991). 

0 + ( 4 S) is the major ion at the F 2 peak. It is produced largely by photoioniza- 
tion of O above 130 km and photoelectron impact ionization at lower altitudes. 
Dissociative ionization of 0 2 is a minor source over the whole ionosphere. 
Charge transfer from N 2 + to O (reaction 249b) is more important at low 
altitudes, where it produces 20-25% of the total source between 100 and 
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150 km. The rate of production of 0 4 ( 4 S) from quenching of 0 + ( 2 D) by N 2 

0 4 ( 2 D) + N 2 > 0 4 ( 4 S) + N 2 (250) 

exceeds that of reaction (249b) above 250 km, if the rate coefficient is close to 
the upper limit of 8 x 10 11 cm 3 s" 1 reported by Johnsen and Biondi (1980). 
Smaller but significant sources of O 4 (on the order of a few percent) include 

N 4 +0 2 ►O 4 +NO (251) 

below about 150 km and 

N 4 +0 >0 4 +N (252) 

above 200 km. No measurements of the rate coefficient for reaction (252) are 
available, but analysis of AE data has yielded values of 2 x 10 _12 cm 3 s _1 
(Constantinides et al , 1979) and 5 x 10“ 13 cm 3 s” 1 (Torr et al , 1979). Bates 
(1989) applied the Landau-Zener approximation to the relevant potential 
curves and computed a value similar to the former. Quenching of 0 4 ( 2 D) by 
electrons 

0 4 ( 2 D) + e ► 0 4 ( 4 S) + e (253) 

produces about one-third of the 0 4 ( 4 S) near 350 km, and quenching by O 
(reaction 175) produces about 5% if the rate coefficient is on the order of 
1 x 10" 1 1 cm 3 s" 1 (Orsini et al, 1977). If reaction (253) is faster, as Oppenhei- 
mer, Dalgarno and Bnnton (1976) have suggested, the source would increase 
proportionately. 

Loss of 0 4 ( 4 S) is mainly by reaction with O a (reaction 243) below 140 km. 
At higher altitudes, reaction with N 2 (reaction 244) is more important. Charge 
transfer to N( 2 D) 

O 4 + N( 2 D) » N 4 + O (254) 

is significant as a loss process for O 4 above about 250 km; it constitutes about 
a third of the chemical loss at 350 km. Of course, above the O 4 peak, near 
300 km, downward transport by diffusion is more important than chemistry in 
determining the O 4 density profile. Altitude profiles of the major chemical 
sources and sinks of O 4 are shown in Figure 1.40. 

0 4 ( 2 D) and 0 4 ( 2 P) are produced mainly by photoionization and electron 
impact ionization of O at all altitudes. Quenching of 0 4 ( 2 P) 

0 4 ( 2 P) + O * 0 4 ( 2 D) + O (255) 

0 4 ( 2 P) + e ► 0 + ( 2 D) + e (256) 

and radiation 

0 4 ( 2 P) ► 0 + ( 2 D) + M7319, 7329 A) (257) 

are the major loss mechanisms for 0 4 ( 2 P) above about 150 km. These processes 
are also secondary sources of 0 + ( 2 D). Quenching of 0 4 ( 2 P) is mostly by O 
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Figure 1.40 Altitude profiles of the major (a) sources and (b) sinks of 0 + for a low 
solar activity model of the terrestrial dayside ionosphere 


below 450 km and by electrons at higher altitudes; loss by radiation exceeds 
that by quenching above 350 km. The reaction of 0 + ( 2 P) with N 2 

0 + ( 2 P) + N 2 > N 2 + + O (258a) 

>0 + ( 2 S) + N 2 (258b) 

>N0 + + 0 (258c) 

is the major loss process for 0 + ( 2 P) below about 150 km. If the measured high 
energy cross sections (Lavolee and Henri, 1989) apply also to thermal energies, 
the major channel is charge transfer (channel 258a). An AE analysis (Oppen- 
heimer et al . , 1976) suggested, however, that at ionospheric temperatures 
quenching (channel 258b) may be more important. Below 350 km, the major 
loss process for O f ( 2 D) is the reaction with N 2 


0 + ( 2 D)+ N 2 - 

— * n 2 + + o 

(259a) 

- 

— - O ’ ( 4 S) + n 2 

(259b) 

— 

— ► NO + + O 

(259c) 
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This reaction proceeds mainly by charge transfer (channel 259a) (Johnsen and 
Biondi, 1980). A minor sink for 0 + ( 2 D) is quenching by O (reaction 175). If 
the rate coefficient for this reaction is less than 1 x 10” 1 1 cm 3 s“ \ it contributes 
a maximum of 3-6% of the total sink from 250-350 km. If the rate coefficient 
is much larger, reaction (175) could be a significant loss process. 

N + is a minor ion in the topside ionosphere of the Earth. Its production is 
mainly through photodissociative and electron-impact dissociative ionization 
of N 2 below 230 km. Direct ionization of N is more important than dissociative 
ionization of N 2 below 230 km. Direct ionization of N is more important than 
dissociative ionization only above about 550 km. Above 230 km, charge transfer 
of 0 + to N( 2 D) (reaction 254) is the dominant source (Constantinides et ai, 
1979). The rate coefficient for reaction (257) was derived from AE data and 
values of 1.3 x 10“ lo cm 3 s _1 were reported by Constantinides et al. (1979) 
and Torr et al. (1979), respectively. Bates (1989) computed the similar value of 
1.4 x 10“ 10 cm 3 s“ l . Charge transfer from 0 4 ( 2 D) to N 

0 4 ( 2 D) 4- N ► N 4 4 O (260) 

may be significant in the altitude range 250-450 km, but the rate coefficient has 
not been measured. A value of 7.5 x 10“ 1 1 cm 3 s” 1 was estimated by Torr et 
al. (1974) from analysis of AE data, and Dalgarno (1979) concluded that the 
data were consistent with an upper limit of 3 x 10“ 10 cm 3 s” 1 . If the rate 
coefficient is 1 x 10“ 10 cm 3 s 1 , reaction (260) contributes about 8% of the 
total source between 350 and 450 km. Dissociative charge transfer of He 4 to N 2 

He + 4 N 2 ► N 4 4 N 4 He (261) 

is a minor source at high altitudes, producing about 10% of the N 4 near 500 km 
(Constantinides et al. , 1979). The products of the reaction of O 24 with N 2 

O 2 4 4 N 2 * products (262) 

are unknown, but even if N 4 is the major product, this reaction provides only 
about 3% of the total source between 150 and 250 km. 

Below about 300-320 km, loss of N 4 is mainly by reaction with 0 2 - The 
yields of channels (244), (241) and (251) are 51, 43 and 6%, respectively 
(Howorka et ai , 1980). At high altitudes, downward diffusion is the dominant 
loss process, but the major chemical loss processes are charge transfer to O 
(reaction 252) and, to a lesser extent, charge transfer to H 

N 4 4 H ► H 4 4 N (263) 

No measurements of the rate coefficients for reaction (263) are available, but a 
value of 3.6 x 10 “ 12 cm 3 s“ 1 was derived from AE data (Constantinides et al. , 
1979) and quantal calculations yield 1.2 x 10” 12 cm 3 s“ 1 at 1000 K (Butler and 
Dalgarno, 1979). 
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Measurements made by the magnetic ion mass spectrometer on the Atmo- 
sphere Explorer satellite revealed 0 2+ number densities ranging from about 
3 cm“ 3 near 150 km to 3CL40 near 400 km (Breig el ai , 1977). The sources of 
0 2+ include photoionization of 0 + , Auger ionization of O by X-rays and 
double photoionization of outer shell electrons. Although measured cross 
sections were not available, Victor and Constantinides (1979) derived plausible 
cross sections for the latter process by analogy to those measured and calculated 
for noble gases. Using these cross sections they showed that double photoioniza- 
tion of valence electrons is the major source of O 2 + in the terrestrial ionosphere. 
The upper and lower limits to the optically thin double photoionization 
frequency for their derived cross sections were about 2 x 10" 9 and 1 x 
10“ 9 s“‘, respectively. Recent measurements (Angel and Samson, 1988) have 
shown that the actual double ionization cross sections are near the lower limit 
values adopted by Victor and Constantinides. For high solar activity, the 
optically thin photoionization frequency implied by the measured cross sections 
is about 1.6 x 10“ 9 cm 3 s l . 

The major loss process for 0 2+ is reaction with N 2 (reaction 262); at low 
altitudes reaction with 0 2 

0 2+ 4- 0 2 ► products (264) 

is also significant. Reaction (264) has been shown to proceed at near gas kinetic 
rates, but the products have not been identified (Johnsen and Biondi, 1978; 
Howorka el ai, 1979). Since the second ionization potential of O is about 35 eV, 
the exothermicities of these reaction are large and therefore, as for reaction 
(262), dissociative channels probably dominate. The charge transfer reaction of 
0 2+ with O 

0 2+ +0 >0 + +0 + (265) 


may be the major loss process at altitudes greater than 350-400 km, but the 
rate coefficient has not been measured in the laboratory. Victor and Con- 
stantinides (1979) derived a crude upper limit of 1 x 10“ 10 cm 3 s” 1 from their 
analysis. Also using AE data, Breig, Torr and Kayser (1982) derived a value of 
6.6 x 10“ 1 1 cm 3 s _1 for a low solar activity optically thin double photoioniza- 
tion frequency of 2.6 x 10“ 9 cm 3 s _l . Analysis Pioneer Venus data yielded a 
rate coefficient for reaction (265) of 1.5 x 10“ 10 cm 3 s“ 1 (Fox and Victor, 1981). 
These models were constructed before experimental values for the double 
photoionization cross sections were available, and should be reconsidered. A 
measurement of the rate coefficient for reaction (265) would be desirable. 

At the high pressures and low temperatures of the terrestrial D region, the 
terminal ions 0 2 + and NO 1 react to form clusters with N 2 , 0 2 and C0 2 , in 
reactions such as (Ferguson et ai , 1979) 


NO + 


+ N, + M 


NO + .N 2 + M 


(266) 
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Then exchange reactions occur, in which the clusters that form are replaced 
with more stable clusters 

NO + .N 2 + H 2 0 > N0 + (H 2 0) 4 N 2 (267) 

Addition of another water molecule 

N0 + (H 2 0) + H 2 0 + M » N0 + (H 2 0) 2 + M (268) 

is followed by a reaction in which a hydrated proton forms 

N0 + (H 2 0) 2 4- H 2 0 ► H + (H 2 0) 2 + HN0 3 (269) 

Proton hydrates H + (H 2 0)„, with n up to four, have been observed by rocket- 
borne mass spectrometers (Kopp and Hermann, 1984). 

Ablation of meteorites deposits metal atoms such as Fe, Al, Ni, Ca and Na 
at altitudes from 80 to 100 km (Zhinden et al . , 1975). Although the fraction of 
Na in meteors is only 0.6% (Mason, 1971), sodium was the first metal atom 
identified in the mesosphere through its intense emission in the nightglow at 
589 nm (Chamberlain, 1961). The sodium layer is located between 86 and 
100 km, with a peak near 93 km of about 3 x 10 3 cm 3 (Plane, 1991). The total 
column density is about 5 x 10 9 cm -2 , but the variability is large, about a 
factor of 10. Metal atoms have low ionization potentials: that for Na is about 
5.1 eV. Metal ions can be formed by photoionization and charge transfer from 
0 2 + and NO + 

0 2 + + Na ► Na + 4- 0 2 (270) 

Na + may then combine with N 2 , forming clusters 

Na + 4 N 2 4 M ► Na\N 2 4 M (271) 

which are replaced with more stable clusters such as Na + .C0 2 and Na + .H 2 0, 
in reactions that are analogous to those of 0 2 + and NO + in the D region 
(Richter and Sechrist, 1979; Thomas et al 1983). These clusters are eventually 
destroyed by dissociative recombination 

Na\X 4 e > Na 4 X (272) 


1.5 SUMMARY 

In this review we have attempted to summarize the ionic reactions that control 
the distribution of ionization in astrophysical and atmospheric environments. 
Despite the large number of reactions that we have described, we may well 
have omitted some that are significant and we have probably included many 
that are not. It is clear the ion chemistries of the different environments have 
elements in common, yet with distinctive differences. Many significant reactions 
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have yet to be studied in laboratory or theoretical studies. Their investigation 
is particularly challenging because of the need often to measure the influence 
of rotational and vibrational populations and to determine the identity and 
state of excitation of the end products of reactions over a wide range of 
temperatures. Important progress has been made, but much remains to be done. 
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Note added in proof 

Recent models of the ionosphere of Triton (J. Lyons, Y. Yung and M. Allen, 
1994, Icarus , in press) suggest that the major ion is C^ produced by photo- 
ionization of C and charge transfer. The most important source of C is 
dissociative recombination of CO + . The electron density profile measured by 
Voyager can be accounted for entirely by solar radiation. 

Although the presence of high mass-28 ion densities in the nightsidc Venus 
ionosphere has been suggested as an indication of the presence of strong 
electron precipitation, recent auroral models indicate that the very high mass-28 
ion densities cannot be reproduced. An inspection of the instrument operating 
voltages, and numerical modeling of the instrument response have suggested 
that offsets in the PV OIMS voltage settings can cause mass-30 ions to be 
collected in the mass-28 setting. Thus the extremely high mass-28 ion densities 
may be an instrumental effect (J. Fox, J. Brannon, J. Grebosky and H. Porter, 
1994, EOS, submitted). 



